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ABSTRACT 


The  equations  of  motion  for  a  ducted  fan  hovering  device 
are  developed  and  programmed  in  a  computer  simulation. 
Experimental  aerodynamic  data  is  integrated  with  the 
computer  model .  A  feedback  control  scheme  for  the  multiple- 
input,  multiple-output  system  is  determined  using  optimal 
control  techniques.  Time  response  results  are  obtained  and 
analyzed. 

As  a  separate  issue,  the  body  of  the  device  is  modelled 
for  electromagnetic  analysis  and  a  basic  antenna  design  is 
determined  for  UHF  transmission. 


Accesiyn  For  ^ 

NT  IS 

CRAil 

D]  1C 

tab  □ 

Un.j'ir.O'jiiccd  □ 

JtiV-tif -C 

.T.ion 

Bv 

DistriL' 

itlO  ;/ 

A 

‘",iy  ..  I 

AVri!  v'Mfl.'Or 

Oist 

h 

i 

1 

j 

3 


TABLE  OF  CONTENTS 


K 


I.  INTRODUCTION . '.  .  11 

A.  WHAT  IS  AROD? .  11 

B.  *  AROD  AND  OPTIMAL  CONTROL . 12 

C.  AROD  AND  RADIO  FREQUENCY  TRANSMISSIONS  ....  13 

D.  SCOPE .  14 

II.  AROD  DYNAMIC  SIMULATION .  16 

A.  THE  GYROSCOPIC  AROD .  16 

B.  EQUATIONS  OF  MOTION .  21 

1.  Angular  Momentum .  22 

2.  Force  Equations .  23 

3.  Kinematic  Equations  of  Motion .  25 

4.  Summary  of  Equations  of  Motion .  27 

C.  APPLIED  FORCES  AND  MOMENTS .  28 

1.  Aerodynamic  Characteristics  .  28 

2.  Forces  and  Momer'  s  for  Control .  33 

3.  Servo  Equations  for  Control  Surfaces  and 

Throttle .  36 

D.  ASSEMBLING  THE  MODEL:  PROGRAMMING .  38 

1.  The  Complete  Block  Diagram .  38 

2.  Program  Flow  Chart .  38 

E.  SIMULATION  RESULTS .  38 

1.  Zero  Input  Response .  41 


4 


<1 


2.  Aileron  Input .  41 

3.  Elevator  and  Rudder  Input .  44 

4.  Throttle  Input .  48 

5.  Response  Without  Aerodynamic  Moments  ...  48 

III.  STABILITY  AND  CONTROL  OF  AROD . .  .  53 

A.  SUMMARY  OF  THE  CONTROL  PROBLEM .  54 

1.  Summary  of  the  Complete  AROD  System  .....  54 

2.  Optimal  Feedback  Control  .  56 

B.  LINEARIZING  THE  AROD  MODEL .  60 

1.  The  Small-Disturbance  Theory  .  . .  61 

2.  Physical  Approximations  of  Forces  and 

Moments . 63 

3.  The  Linear  Model  Assembled .  65 

C. .  FORMULATING  THE  CONTROL  PROBLEM .  66 

1.  State  space  Representation  .  66 


2.  Stability,  Observability,  and 
Controllability  . 


3.  Discretizing  the  System .  75 

D.  THE  OPTIMAL  REGULATOR . 79 

1.  Optimal  Control .  80 

2.  Selecting  the  Performance  Measure  .  82 

4.  Iterative  Design  Process  .  85 

5.  Linear  System  Optimal  Regulator  Response  .  89 

E.  THE  OPTIMAL  TRACKER .  94 

F.  RESULTS  WITH  THE  NONLINEAR  SYSTEM .  95 


1.  Nonlinear  Model  without  Aerodynamic  Moments  99 


5 


2.  Nonlinear  Modal  Complete 

3 .  Conclusion  . 


103 


.  109 


IV.  RADIO  FREQUENCY  ANTENNA  DESIGN  .  110 

A.  REQUIREMENTS  FOR  A  RF  ANTENNA . 110 

B.  DESIGN  SELECTION . .  .  112 

C.  THEORETICAL  DEVELOPMENT . 114 

*  1.  Dipole  Array  . . .  .  114 

2.  Loop  Antennas . 116 

3.  Ground  Plane  Effects  .  118 

D.  BODY  MODELLING  AND  COMPUTER  AIDED  DESIGN  RESULTS  121 

1.  Results  with  Dipole  Arrays . 122 

2.  Results  with  Loop  Antennas . 130 

E.  ANALYSIS  &  DESIGN  CONCLUSIONS . 138 

V.  CONCLUSIONS . 141 

LIST  OF  REFERENCES . 143 

APPENDIX  A  EXPERIMENTAL  DATA . 145' 

APPENDIX  B  SIMULATIONS  SOURCE  CODE . 149 

APPENDIX  C  CONTROLLER  AND  LINEAR  MODEL  SIMULATION  .  .  158 

APPENDIX  D  NUMERICAL  ELECTROMAGNETICS  CODE  .  161 

1 


i  i  . . . .  . .  mm  mmmm  miiTM 


TABLE  QT  SYMBOLS 


SxmOSl  Description 

Coordinate  Definitions 

X,  Y ,  Z  *'  Body-Fixed  Coordinate  Axes 

\  ,j ,k  Unit  vectors  in  X,Y,Z,  respectively 

X ' ,  Y ' ,  Z '  Earth-Fixed  Coordinate  Axes 


Dynamic  States 


tf,v,w 

U' ,V' ,W' 
FrQ.R 


u 

wr 


Body-Fixed  Velocities  along  X,Y,Z 
Earth-Fixed  Velocities 
Body-Fixed  Roll,  Pitch,  and  Yaw 
Earth-Fixed  Euler  Angles 
Total  Body  Angular  Velocity 
Rctor  Angular  Velocity 


Cgna&anSa 

lx  * 

Ir 

m 

g 

lacgfla 

Fx  ,  Fy  ,  F, 
x  »  y  »  F|  j 
X  »  Ft  y  •  z 
F,  ,Ft  ,  Fd 


Body  Moments  of  Inertia  (Body  Axes) 
Rotor  Moment  of  Inertia  (Rotor)  Axis 
Total  Mass  of  Vehicle 
Gravitational  Constant  (32.174) 

Total  Translational  Forces  (Body) 
Aerodynamic  Translational  Forces 
Thrust  and  Propeller  Induced  Forces 
Side-Slip,  Lift,  and  Drag  Forces 


7 


Units 


ft/ sec 
ft/ jec 
rad/sec 
radians 
rad/sec 
rad/ sac 

ft2  -lb  (m) 
ft*  -lb  (m) 
lb(m) 
ft/ sec2 

lbs 

lbs 

lbs 

lbs 


Monanta 
L,M,N 
L.  ,M.  ,N, 
L(  ,Mt  ,Nt 
Rr  , Pp 


Total  Angular  Monants  about  X,Y,Z  ft-lba 


Aarodynanlc  Angular  Monants 


ft-lba 


Thrust  and  Propallar  Znducad  Monants  ft-lbs 
Rol 1 , Pitch , Yaw  ( ralativa  to  Vt  0 1 )  ft-lbs 


Angular  Nonantun  of  Body 


ua  ,  ,  ur ,  ut  Input  Signals  to  Sarvo 

s,  ,St  ,sr  ,St  Output  fron  Sarvo 


ft-lbs/sac 


volts 


radians 


Engina  Thrust  Fores 


Vectors 


x(t)/x(k)  Continuous/ Discrata  Tina  Systan  States 
y (t)/y (k)  Continuous/ Discrata  Tina  Output  Statas 
rc  (t)/re  (k)  Continuous/Discrata  Tina  Rafaranca  Input. 
Uc  (tj/u,.  (k)  Continuous/Discrata  Tina  Control  Input 


Matrices 


Continuous  Tina  Systan  Plant 

Continuous  Tine  Control  Distribution 

Discrete  Tine  Systan  Plant 

Discrata  Tine  Control  Distribution 

Output  Measurement 

Feedback  (continuous  and  discrata) 

Observability 

Controllability 

Final  State  Weighting 


Transient  State  Weighting 


8 


R  Control  Coat  Weighting 

J  Cost  Function  (Scalar) 

Sampling  Quantitiaa 

AT  Sampling  Interval  (Period) 

f.  Sampling  Frequency 

-n  System  Natural  Frequency 

Sttlfilftl  St&tM 

Body-Fixed  Pitch  and  Yaw  Angles 


9 


seconds 
herts  (Hz) 
hertz  (Hz) 

radians 


JLuai-^a.'  lauai  wu 


iiaiU 


Especially  to  Shirley,  my  wife,  whose  encouragement  and 
support  was,  is,  and  will  always  be  sincerely  appreciated. 


10 


I.  INTRODUCTION 


A.  WHAT  IS  AROD? 

Tha  U.  s.  Marin*  Corps  haa  aatablishad  tha  Ground/Air 
Tala  Robotic  Systems  (GATERS)  program  to  develop,  tast  and 
field  robotic  devices  fcr  diract  usa  by  U.  S.  Marine 
infantryman.  These  robotic  systems  will  be  operated  under 
the  concept  of  c ale  robotics  which  implies  that  a  man, 
stationed  remotely  from  tha  vehicle,  will  control  the  device 
as  if  ha  is  actually  on  board. 

Tha  Airborne  Remotely  Operated  Device  (AROD)  is  the 

vehicle  being  developed  through  the  GATERS  program  to  meat 

the  Marin*  Corps'  needs  for  an  airborne,  short-range,  direct 

support  remotely  piloted  vehicle  (RPV) .  The  GATERS  program 

office_[l]  states  that,  from  a  tactical  perspective: 

The  AROD  is  ...  designed  to  allow  the  front  line  commander 
to  see  "over  the  next  hill",  out  to  a  distance  of  two 
kilometers  and  "around  the  next  corner"  in  an  urban 
environment . 

This  tactical  requirement  has  led  to  a  more  technical 
statement  of  what  AROD  will  be.  The  GATERS  program  has 
stated  [ 2 ] : 

The  AROD  will  be  a  ducted-fan,  hovering  device  capable  of 
a  forward  translation  speed  of  30  miles  per  hour.  The 
vehicle  will  carry  a  fibre  optic  data  link  and  on  board 
cameras  to  support  assigned  missions  out  to  a  distance  of 
5  kilometers.  Normal  operating  radius  is  2  kilometers. 

More  detailed  requirements  are  listed  in  Appendix  A. 
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B.  APOD  AND  OPTIMAL  CONTROL 


The  AROD  is  to  be  operated  by  a  Marine,  untrained  and 
unskilled  in  aerodynamic  flight.  Therefore,  .  AROD  must  be 
inherently  stable  and  responsive  to  simple  controls  such  as 
"go  up/down",  "rotate  right/left" ,  "go  forward",  "go 
sideways" .  These  instructions  will  be  commanded  from  a 
joystick,  through  the  control  uplink  to  an  on-board  control 
system.  The  control  system  will  interpret  the  command  and 
perform  the  necessary  tasks  to  maneuver  the  AROD.  The 
design  of  this  controller  is  the  major  topic  of  this  work. 
Figure  1.1  is  a  schematic  of  the  prototype  AROD. 


Figure  1.1  Airborne  Remotely  Operated  Device  (AROD) 
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Chapter  II  discusses  the  AROD  system  and  the  equations 
which  govern  its  motion.  The  equations  are  complex  and  the 
system  has  multiple  inputs  and  multiple  outputs.  One 
approach  to  controller  design  for  such  a  system  is  to 
simplify  the  equations  describing  the  system  as  much  as 
possible  while  retaining  as  much  of  the  system  characteris¬ 
tics  as  possible.  This  is  the  method  developed  herein.  ~ 

The  complex  multiple-input,  multiple-output  AROD  system 
is  appropriately  addressed  using  optimal  control  techniques 
with  the  aid  of  a  digital  computer  [3].  Optimal  control, 
therefore,  is  the  approach  taken  and  developed  in  this  work 
to  solve  the  problem  of  stability  and  control  of  the  AROD. 

C.  AROD  AND  RADIO  FREQUENCY  TRANSMISSIONS 

AROD,  as  previously  described,  utilizes  a  fiber  optic 
tether  to  pass  control  information  and  a  video  signal  from 
on  board  observation  devices.  Therefore,  the  need  for  radio 
frequency  (RF)  transmissions  to  or  from  the  vehicle  is  not 
immediately  obvious.  In  fact,  under  ideal  tactical  condi¬ 
tions,  the  RF  signature,  susceptible  to  hostile  electronic 
countermeasures,  would  be  a  most  undesirable  characteristic 
for  a  flying  vehicle  to  have.  Nevertheless,  the  fiber  optic 
tether  will  be  subject  to  damage  and  a  back-up  system  must 
be  considered  for  vehicle  control,  and  for  sending  back  the 
video  signal  as  well.  Most  importantly,  training  must  be 
conducted  prior  to  any  conflict  or  situation  in  which  AROD 
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is  to  be  used.  The  fiber  optic  cable,  not  immediately 
reusable  by  the  lower  echelon  ground  forces,  is  expensive 
and,  again,  susceptible  to  damage.  Therefore,  it  is 
desirable  to  provide  a  low  cost  method  for  training  opera¬ 
tors  and  tactical  commanders  in  the  use  of  the  AROD;  hence 
the  need  for  an  alternative  to  a  fiber  optic  tether. 

A  back-up  RF  communications  link  to  provide'  operator 
commands  to  the  AROD  has  already  been  designed  and  imple¬ 
mented.  However,  this  is  an  uplink  only  and  provides  no 
means  for  the  on  board  cameras  to  return  their  wide 
bandwidth  video  signal  to  the  ground  station.  Some 
consideration  of  this  need  has  been  made  in  the  GATERS 
program  and  a  transmission  frequency  in  the  800  to  900 
megahertz  range  has  been  tentatively  approved. [2]  It  has 
also  been  recognized  that  a  video  signal  via  RF  might  be 
adequate  even  though  it  was  somewhat  degraded  relative  to 
the  high  bandwidth  video  signals  possible  with  fiber  optics. 
With  these  assumptions  in  mind.  Chapter  IV  deals  with  the 
design  of  an  on  board  antenna  suitable  for  transmitting  a 
video  signal  from  AROD  to  the  ground  station. 

D.  SCOPE 

In  summary-,  this  work  will  address  two  different 
problems  associated  with  AROD: 

(1)  Design  of  a  control  system  for  stable  flight. 

(2)  Design  of  a  radio  frequency  antenna  for  the  AROD. 
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1.  Controller  Design 

System  design  of  the  vehicle  has  been  previously 
undertaken  by  AROD  project  engineers.  The  prototype  AROD 
has  been  built  and  pertinent  information  is  listed  in 
Appendix  A.  In  Chapter  II ,  the  model  describing  AROD  is 
developed.  Design  of  the  proposed  control  system  is  under¬ 
taken  in  Chapter  III. 

2,  Antenna  Design 

Chapter  IV  deals  exclusively  with  design  and 
analysis  of  an  RF  antenna  proposed  for  AROD.  The  design 
approach  is  to  analyze  the  requirements  and  determine  a 
basic  design  which  can  be  refined  gradually  until  an 
acceptable  solution  is  found.  The  computer  program 
Numerical  Electromagnetics  Code  (NEC) ,  developed  by  the 
Naval  Ocean  Systems  Center  is  used  extensively  to  aid  in  the 
design.  No  attempt  is  made  to  recommend  specific  hardware 
components  in  order  to  implement  the  antennas. 
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II.  AROD  DYNAMIC  SIMULATION 


The  theory  behind  the  equations  governing  aerodynamic 
flight  is  well  developed.  Procedures  exist  which  make  the 
modelling  of  complicated,  multiple  engine,  highly  maneuver- 
able  aircraft  almost  a  step  by  step  task.  Unfortunately,  the 
dynamics  which  characterize  AROD  differ  considerably  from 
the  standard  flight  model  at  a  very  basic  level.  Like  all 
flying  vehicles,  AROD  has  six  degrees  of  freedom  of  motion: 
three  linear  displacements  and  three  angular  rotations,  but 
here  the  similarities  end.  This  ducted  fan  device  has  no 
lift  prod  icing  surfaces  other  than  the  propeller  which 
accounts  for  100  percent  of  its  lift  and  100  percent  of  its 
complexity.  This  chapter  traces  the  development  of  a 
suitable  dynamic  model  of  the  AROD,  concentrating  first  on 
the  significant  gyroscopic  contribution  of  its  propeller, 
and  then  on  the  complete  vehicle.  The  goal  is  to  provide  a 
computer  simulation  with  which  an  experimental  control 
algorithm  may  be  integrated  for  the  purpose  of  analysis  and 
improvement.  [4:pp.2-22] 

A.  THE  GYROSCOPIC  AROD 

AROD  is  a  gyroscope.  The  single  propeller  rotates  about 
the  longitudinal  vehicle  axis  to  produce  a  downwash  or  jet 
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of  air  through  the  duct  which  makes  up  the  AROD  body.  This 
jet  is  directed  by  pitching  the  body  forward,  backward,  or 
side  to  side  so  that  the  vehicle  moves  away  from  the  jet  in 
the  direction  desired.  To  appreciate  tho  control  effort 
required  pitch  the  vehicle  in  a  specific  direction 

without  causing  undesired  effects  (coupling)  in  the  other 
five  degrees  of  freedom,  the  equations  governing  a  gyroscope 
need  to  be  understood. 

Figure  2.1  pictures  a  gyroscope  which  rotates  in  the 
same  plane  as  the  AROb  propeller  relative  to  the  coordinate 
system  described.  This  axis  orientation  will  also  serve  to 
define  symbols  for  later  use.  For  illustration,  the  body- 
fixed  and  earth-fixed  axes  do  not  coincide  in  Figure  2.1. 
X, Y, Z  and  P,Q,R  are  the  body-fixed  descriptive  geometry  for 
linear  displacement  and  angular  rotation  rate.  P,Q,  and  R 
will  also  be  referred  to  as  roll,  pitch,  and  yaw.  X',Y',Z' 
and  are  the  corresponding  earth-fixed  displacements 

and  rotations.  X'  can  be  thought  of  as  altitude.  p,6,  and 
$  are  the  Euler  angles.  Velocities  will  be  denoted  by  U,V, 
and  W  for  rate  of  displacement  in  the  X,Y,  and  Z  directions, 
and  U' , V'  ,W  for  X*  ,Y*  ,Z* .  Pictured  also  for  use  in  later 
discussion  are  L,M,N  (rolling,  pitching,  and  yawing  moments) 
and  Fx,Fy,Fz  (applied  force  components). 

Angular  momentum,  h,  due  to  the  propeller  is  defined  as 

hr  *  Ir wr  *  i  hr x  +  j  hP y  +  k hP z  (2.1) 
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The  rotor  lies  in  the  YZ  plane  of  the  body  coordinates. 
Since  the  coordinate  axes  have  been  chosen  so  that  the  rotor 
spins  symmetrically  about  the  X-axis,  hr  is  directed  only 
along  X  and  hPy«  hP4-  0.  Equation  2.1  becomes 

hr  -  Ir«P  -  fhri  (2.2) 

At  this  point,  a  simple  calculation  will  illustrate  the 
significant  influence  the  rotor  has  over  the  AROD.  Refer¬ 
ence  4  states  that  in  developing  the  equations  of  motion  for 
aircraft  with  rotors  which  exert  gyroscopic  moments  on  the 
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body  "more  of  ton  than  not,  such  gyroscopic  moments  turn  out 
to  be  negligible."  In  the  AROD  case,  however,  Ir  -  .015 
ft2 -lb (mass)  and  hover  speed  of  the  propeller  is  wP  -  7200 
rpm  -  753.98  rad/ sec.  From  Eq.  2.2,  hr-  11.3  ft2 -lb (m) /sec. 
Compared  with  a  nominal  total  mass  of  2.64  lb(m)  (85  lb 
weight) ,  it  is  clear  that  the  angular  momentum  imparted  by 
the  rotor  is  significant  and  that  gyroscopic  effects  'will 
play  a  large  part  in  modelling  the  dynamic  behavior  of  AROD. 

Complications  arising  from  these  gyro  effects  can  also 
be  illustrated.  If  the  AROD  body  and  rotor  are  considered  a 
complete  gyroscopic  system,  then  the  angular  momentum  of  the 
body  must  be  included  in  the  formulation.  Angular  momentum 
due  to  a  rotating  body  is  noted  ash*,  —  i  hb  x  +jhby  +khb2. 
With  the  orientation  of  Figure  2.1,  Eq.  2.1  becomes 

h  -  \  (h*,*  +  hTK)  +  jhfcy  +  (2.3) 

-  i  (PI*  +  lr«r)  +  jQIy  +  kRl* 

where  h  is  the  combined  angular  momentum  of  the  rotor 
and  body. 

Angular  moments  about  the  axes  of  a  gyro  are  discussed  in 
the  next  section  and  defined  by  Euler  [5:p.  93]  as 

L  -  h*  -  hyR  +  h*Q  (2.4) 

M  -  hy  -  h*P  +  h*R 
N  -  h,  -  h*Q  +  hyP 

A  simulation  of  this  simple  gyro  model  is  given  in 
Appendix  B.  Results  are  shown  in  Figure  2.2  for  the  case  of 


mu 
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a  disturbance  about  the  Y  axis  (pitch) .  Hots  that  the 
actual  moments  of  inertia,  listed  in  Appendix  A,  wore  used 


Figure  2.2  Gyroscopic  Coupling  Effects 

in  this  first  step  computer  simulation.  The  time  response 
shown  here  is  for  the  simple  case  with  no  rotor  accelera¬ 
tion  term  («r  *0)  and  is  the  result  of  a  one  secon  step 
for  an  angular  moment,  M,  beginning  at  time  -  1  second.  The 
amplitude  is  unimportant  and  has  no  effect  on  the  time 
response  except  to  scale  the  result.  The  purpose  of  the 
input  was  to  pitch  over  the  coordinate  system  by  applying  a 
torque  about  the  Y  axis  (Q)  .  Ideally,  a  stable,  well 
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behaved  device  would  simply  pitch  over  and  maintain  its  naw 
attituda.  This  is  not  tha  casa  for  AROD  and  Figure  2.2 
illustratas  tha  pracassion  and  nutation  duo  to  a  savara 
coupling  of  tha  pitch (Q)  and  yaw(R)  axes.  Thasa  gyroscopic 
affects,  than,  ara  tha  reason  a  "simple"  davica  such  as  AROD 
prasants  such  a  complex  controls  problem  to  tha  designer. 


B.  EQUATIONS  OF  MOTION 

Thi  basic  aquations  of  motion  ara  developed  for  a 
typical  aircraft  in  References  4  and  5.  A  combination  of 
the  two  approaches  is  taken  here  to  arrive  at  tha  sat  of 
equations  describing  AROD. 

The  Euler  aquations  of  motion  for  a  free  flying  rigid 
vehicle  with  six  degrees  of  freedom  are  derived  from 
Newton's  second  law.  For  angular  moments  applied  to  the 
coordinate  axes,  we  have 


L 

N 

N 


,4ll 

’dt 

Jdt 

*dt 


(2-5) 


where  h»  i  hx  +  j hy  +  khz  and  is  the  vehicle's  total 
angular  momentum. 


For  linear  or  translational  forces  applied  along  the  axes. 


F* 


<5 


1  Vt 


F, 


(2.6) 
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where  v  ■  i v  +  j v  +  kv  and  is  tha  total  linaar  transla¬ 
tional  valocity  of  tha  cantar  of  mass(m)  of  AROD. 

Tha  dynamic  charactariatics  of  AROD  ara  raprasantad  by 
Equations  2.5  and  2.6.  Othar  ralationships  ara  required  to 
solva  for  aach  unknown  and  for  aach  physical  forca  and 
momant  acting  on  tha  body.  Tha  davalopaant  of  a  coaplata 
modal  is  undartakan  now  for  tha  apacific  casa  of  AROD  using 
actual  mcasuramants  and  experimental  data  from  tha  proto- 
typa. 

1.  ftnmiac  Moaantua 

Equation  2.4  is  axpandad  to  bacoma 


whara 


h*-  i  (IXP  +  Irwr)  +  j  (Ix PR  ►  IPwPR) 

-  k  (I*  QP  +  Ir  wr  Q) 

fay-  -1  IyRQ  +  j  IyQ  +  kIyPQ 
h,“  J  IXQR  -  j  Ix PR  +  k  ItR 

&  *  ft  +  "  x  h 

h  -  i  (IXP  +  IrwP)  +  jlyQ  +  kItR 
w  -  i  P  +  J  Q  +  kR 


(2.7) 


Combining  Equations  2.4  and  2.7,  tha  AROD  angular  momant 
equations  are : 


L  -  IXP  +  (Is  -  Iy)QR  +  IPwP  (2.8) 

M  -  IyQ  +  (Ix  -  IX)PR  +  RIrwr 

N  -  IXR  +  (Iy  -  IX)PQ  -  QIr wP 


The  rolling,  pitching  and  yawing  moments  (L,  M,  and  N)  may 
be  induced  by  the  aerodynamics  of  the  body  and  environment 
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or  by  control  inputs  through  ths  control  vanss  and  ars 
rafarrad  to  as  tha  "appliad  moments".  Thay  ara  discussad  in 
a  later  sectirn.  Equation  2.8  is  next  rearranged  so  that 
tha  unknown  variables  (dynamic  staves  of  interest)  nay  be 
determined.  Roll,  pitch,  and  yaw  rates  are: 

P  -  ^C<ly  -  I«)Q*  -  Ir«r  -  L]  .  (*^> 

Q  *  ^[(i*  "  -  it] 

R  -  ISC(I*  -  Iy)PQ  +  QIr«r  *  »] 

Already  obvious  is  the  mathematical  explanation  for 
the  coupling  that  was  observed  in  the  last  section.  Every 
combination  of  pitch,  roll,  and  yaw  are  present  in  Equation 
2.9.  Pitch,  roll,  and  yaw  are  each  "coupled"  to  the  others. 
Additionally,  coupling  of  the  rotor  angular  acceleration 
with  roll  will  cause  the  throttle  controlling  the  propeller 
speed  to  influence  these  states  and  further  complicate  the 
system. 

2.  gores.  Equation! 

The  expressions  describing  the  translational  forces 
are  developed  from  the  body  accelerations  which  are 

V„-  IU  +  jRU  -  kQU  (2.10) 

vy-  -i  RV  +  J  V  +  k  PV 
vt-  i  QW  -  j  FW  +  k  W 

where  v»  iv„  +  jvy  +  kv,  and  is  the  total  body  velocity 
vector. 
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Equations  2.6  and  2.10  ara  oosbinad  and  raazxanqad  to  obtain 
tha  differential  aquations  containing  forca  tarns, 

6  «  (RV  -  QW)  ♦  (2.11) 

v  -  (pw  -  r a)  +  Jy 

W  «  (QO  -  FV)  +  **  _ 

whara  n  is  tha  nass  of  tha  vahicla. 

Equations  2.11  contain  nonlinear  tarns  as  did 
Equation  2.9.  Thasa  nonlinaaritias  ara  in  tha  form  of 
products  of  status.  Tha  forca  coaponants  F„ , Fy , F,  consist 
of  lift  from  tha  propaller  (Ft ) ,  other  aerodynamic  influ¬ 
ences  (Fa ) ,  and  gravity  (g) .  Ft  and  Fa  depend  largely  on 
experimental  data  in  this  development  and  are  taken  up 
later. 

-  -  Gravity  is  a  constant  acceleration  causing  a 
constant  force  in  tha  -X'  direction  (earth- fixed) .  Note 
that  small  changes  in  weight  due  to  fuel  loss  are  neglected. 
As  the  body  rotates,  the  body-fixed  axes  change  relative  to 
the  earth-fixed  axes.  Therefore,  a  coordinate  translation 
is  needed  to  relate  gravity  to  each  of  the  X,Y,Z  axes 
(9x»9y#9z) •  The  gravity  components  are  given  by 

» •  g-  »g«  +  i9y  +  kg*  (2.12) 

Any  translation  from  earth-fixed  to  body-fixed  axes  can  be 
described  by  a  rotation  first  about  the  Euler  angle  p,  then 


9,  and  lastly  Figure  2.3  illuatratas  thaaa  thraa 

rotations  and  tha  rasulting  ralationships  batwaan  coordinata 
systems.  From  Flgura  2.3  and  Bq.  2.12,  tlia  gravity  com- 
ponants  are 

gE  ■  g(cosy»  cose)  (2.13) 

gy  -  -g(sin*  cose) 
g,  -  g(«ine) 

Tha  first-order  diffarantial  aquation  for  translational 
valocitias  can  now  ba  writtan  as 

U  -  (RV  -  QW)  +  gE  +  (2.14) 

V  -  (PU  -  RO)  +  gy  + 

if  -  (QU  -  PV)  +  gE  +  ^ 

We  now  hava  aquations  dascribing  tha  six  dagraas  of 
freedom  of  AROD  and  tha  foundation  for  tha  nodal  has  baan 
laid. 

3.  Klneaftfcic  Equationa  at  Mftfclfln 

Tha  coordinata  translation  through  tha  Eular  anglas 
rasulted  in  an  expression  for  relating  gravity  to  the 
rotating  body.  However,  the  price  paid  for  this  necessary 
convenience  is  the  addition  of  three  more  variables  (y>,6,^) 
to  the  equations.  Therefore,  three  more  equations  are 
needed  and,  if  possible,  the  three  variables  should  be 

expressed  in  terms  of  the  existing  force  or  moment  vari¬ 
ables.  Fortunately  a  relationship  can  be  found  in  the 


-sin*>  cos© 


equation  describing  the  body  angular  rotations.  Angular 
movement  of  AROD  has  been  expressed  in  terms  of  the  body- 
fixed  angular  rates  P,Q,R.  It  can  also  be  described  in 
terms  of  the  Euler  rates  <p,  6,  and  ^  as 

u>  •  I  P  +  jQ  +  kR  “  \  '  <p  +  j'0  +  k'0  (2.15) 


“A  coordinate  transformation  similar  to-  the'  one 
illustrated  in  Figure  2.3  is  possible  to  relate  pitch,  roll, 
and  yaw  with  the  Euler  angle  rates.  However,  some  simpli¬ 
fications  can  be  made  with  each  rotation  since  are,  in 
fact,  the  variables  of  interest.  The  resulting  transforma¬ 
tion  is  called  the  kinematic  eguations  of  motion: 


1  sinp  co&ptanO 

‘p“ 

0  cosp  -sinp 

Q 

0  sing  gftSig 

sine  cose 

XV 

•  « 

.  m 

4.  Suwmar-Y  of  Equations  of  Motion 


(2.16) 


Nine  equations  and  nine  unknowns  have  been  developed 
that  completely  describe  the  dynamic  behavior  of  AROD. 
P,Q,R,  and  U,V,W  represent  the  six  degrees  of  freedom  and 
relate  the  constant  force  exerted  by  gravity  to  the 
body-fixed  coordinates.  Through  the  Euler  angles,  the 
behavior  of  AROD  can  also  be  observed  relative  to  the  Earth. 
Specifically,  the  translational  velocities  for  the  earth- 
fixed  system  U,,V,,W'  can  be  determined  from  body- fixed 
velocities  U,V,W  and  p,e,y>  using  the  coordinate  transforma- 
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tion  of  Equation  2.13.  The  applied  forces  and  moments  are 
made  up  of  external  forces  and  commanded  inputs,  and  are 
represented  by  Fx  ,  Fy ,  Fz  ,  and  L,  M,  N.  AROD  can  be 
pictured  as  the  block  diagram  in  Figure  2.4.  The  external 
influences  are  taken  up  next. 

C.  APPLIED  FORCES  AND  MOMENTS 

The  applied  forces  and  moments  are  a  result  of  the 
commanded  inputs  and  the  movement  of  AROD  through  the  air. 
The  commanded  inputs  can  change  the  rotor  speed  causing  a 
change  in  the  lift  generated  and  a  torque  on  the  body. 
Commands  can  also  come  in  the  form  of  control  surface 
displacements  which  steer  the  AROD  and  cause  the  body  to 
rotate  about  its  axes.  The  movement  of  AROD  through  the  air 
will  create  drag.  Drag  forces  can  also  apply  both  transla¬ 
tional  forces  and  angular  moments  on  the  body  and  are  a 
result  of  AROD's  aerodynamic  characteristics. 

i.  Asradynanis  fiharasfcgciafcisg 

Two  methods  exist  for  determining  the  aerodynamic 
forces  and  moments  on  a  vehicle  [4sp.lll]:  (1)  by  experiment 
or  empirical  analysis  and  (2)  by  theoretical  development. 
The  great  advantage  of  (1)  is  accuracy  and  the  advantage  of 
(2)  is  low  cost.  The  fact  that  AROD  is  a  ducted  fan  causes 
a  complex  interaction  of  the  rotor  with  the  body,  and  the 
downwash  with  the  control  surfaces.  Hence,  the  ducted  fan 
compounds  the  accuracy  problem  and  empirical  means  for 
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Figure  2.4  Block  Diagram  of  AROD  Equations  of  Motion 
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determining  the  forces  and  moments  is  the  most  suitable 
approach.  Sandia  National  Laboratories  has  built  a  proto¬ 
type  AROD  and  has  collected  data  based  on  wind  tunnel  tests 
of  the  AROD  [6].  This  data  consists  of  tabular  results 
describing  the  aerodynamic  drag  coefficients  and  downwash 

swirl  effect;  stability  derivatives;  and  physical  measure- 

*  .  _  — 

ments  of  constants  such  as  weight,  moments  of  inertia,  and 

servo  gains.  This  data  is  listed  in  Appendix  A.  The 
tabular  data  forms  a  basis  from  which  applied  forces  and 
moments  may  be  determined  once  the  following  flight  informa¬ 
tion  is  known: 

(1)  Weight  ratio  of  actual  AROD  to  85  lb  prototype.  The 
data  was  originally  obtained  with  an  85  lb  prototype 
and  must  be  scaled  for  other  weights  such  as  the 
model  discussed  here  which  is  76.5  lb. 

(2)  Propeller  angular  speed.  Propeller  speed  is 
determined  by  setting  a  hover  speed  (this  model  7200 
rpm)  and  calculating  changes  from  the  throttle 
setting  of  the  engine. 

(3)  Propeller  angular  acceleration.  Acceleration  is 
based  on  the  difference  between  speed  calculated 
from  the  throttle  and  speed  at  hover.  This  is  an 
experimental  approximation  arrived  at  by  AROD 
project  engineers. 

(4)  Translational  speed  U,V,W.  Aerodynamic  forces  on 

the  body  are  based  on  the  speed  and  direction  of 
AROD  (angle  of  attack) . 

The  forces  and  moments  which  result  are  of  two 
types:  aerodynamic  and  thrust.  The  forces  and  moments  due 
to  the  aerodynamic  data  are  referred  to  as  Fax  ,  Fay,  Faz, 
and  La  ,  Ma ,  Na  .  Those  due  to  thrust  and  propeller  con¬ 
siderations  are  Ft  x  ,  Ft  y  ,  Ft  z  ,  and  L*  ,  Mt  ,  Nt  . 
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No  attempt  is  made  in  this  work  to  verify  the  data. 
The  relationships  leading  to  the  forces  and  moments  dis¬ 
cussed  in  this  section  are  taken  verbatim  from  the  informa¬ 
tion  supplied  by  AROD  project  engineers. 

a.  Angle  of  Attack 

The  tabular  data  (aerodynamic  drag  coefficients 
and  downwash  swirl)  result  in  forces  and  moments  relative  to 
the  total  velocity  vector,  Vtot,  of  the  body.  A  relation¬ 
ship  is  needed  to  transform  these  forces  and  moments  to 
body-fixed  axes  forces  and  moments.  Angles  of  attack  (a,/?) 
are  the  standard  relationship  for  this  purpose.  Angle  of 
attack  (a)  is  normally  defined  as  the  angle  in  the  XZ  plane 
between  Vtot  and  the  X  axis.  Angle  of  attack  (0)  is  the 
angle  in  the  XY  plane  between  VtQt  and  the  X  axis  (sometimes 
referred  to  as  the  side  slip  angle). [4]  The  angle  of  attack 
used  in  this  model  is  defined  and  illustrated  in  Figure  2.5. 
The  resulting  relationships  are  given  by 

Bina  -  r® —  (2.17) 

vtot 

sin/?  -  -j1— 

vtot 

where  Vt  0  t  -  ,/  U2+  V2+  W2 

b.  Aerodynamic  Forces 

The  forces  associated  with  the  tabulated  data 
are  lift  (Ft  ) ,  drag  (Fd),  and  side  force  (Fs  )  .  They  are 
also  depicted  in  Figure  2.5.  The  transformation  from  lift, 
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Figure  2.5  Angle  of  Attack  Illustration 


drag,  and  side  force  is  accomplished  in  like  manner  to  the 
transformation  in  Figure  2.3  and  results  in 


"  sinacosl  -cosacosl  ain/T 

v 

F.y 

- 

-sinasinl  cosasinl  cos fi 

Fd 

*.«. 

-cosct  -sina  0 

F.. 

c.  Aerodynamic  Moments 

In  similar  fashion  to  the  aerodynamic  forces, 
the  moments  applied  to  the  body  axes  as  a  result  of  AROD's 
movement  through  air  can  be  derived.  These  moments  are 
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rtf erred  to  as  the  aerodynamic  angular  moments  of  roll, 
pitch,  and  yaw  (L,  ,  M,  ,  N,  )  and  are  given  as 


"  sinncos 0 

-cosacoBl 

sin f 

V 

(2.19) 
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-sinasinl 

cosasinl 

COB0 

Rr 

». . 

-cosa 

-sina 

0 

.PP 

*  -  ■— 

where 

Yy  is  the  yaw  moment  relative  to  Vtot  , 

Rr  is  the  roll  moment  relative  to  vt  0 1 , 

Pp  is  the  pitch  moment  relative  to  Vt  0 1 . 

2.  Forces  and  Moments  for  Control 

The  forces  and  moments  previously  discussed  are 
based  on  observation  of  an  AROD  prototype  and  are  external 
to  tha  system  being  designed.  A  second  category  is  those 
forces  and  moments  which  are  a  direct  result  of  the  com¬ 
manded  inputs.  As  mentioned  earlier,  the  inputs  control  (l) 
the  rotor  speed  and  (2)  the  displacement  of  the  control 
surfaces . 

a.  Moment  Due  to  Ducted  Fan  Effects 

The  AROD  in  hover  (constant  rotor  speed)  closely 
resembles  a  gyroscope.  A  gyroscope  imparts  no  torque  on  its 
axi3  If  it  spins  with  a  constant  angular  rotation.  If  the 
rotw  accelerates  (speeds  up  or  slows  down)  a  torque  is 
applied  to  the  axis.  This  torque  is  accounted  for  in 
Equation  2.8.  However,  AROD  is  also  a  ducted  fan  and  the 
drag  between  the  rotor  tip  and  the  inside  body  wall  creates 
a  moment  about  the  X  axis  (P,  roll) .  The  project  engineers 
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for  AROD  have  determined  an  approximation  for  this  moment 
based  on  experiment  which  is  given  as 

I^r  «  .0729  (Awr2)  (ft-lbs)  (2.20) 

where  A wP  is  the  difference  between  the  rotor  speed  and 
hover,  Subscript  Tr  refers  to  a  dependence  on  the 
throttle,  and  .0729  is  a  constant  due  to  the  duct  and  is 
referred  to  as  in  Appendix  A.  _ 

b.  Moments  Due  to  Control  Surface  Displacement 

Command  inputs  also  include  the  ability  to 
displace  the  four  control  vanes  into  the  downwash  from  the 
duct.  As  these  surfaces  are  displaced,  they  impart  moments 
about  the  body  axes.  The  vanes  are  arranged  symmetrically 
as  in  Figure  2.6.  The  vanes  are  displaced  by  a  servo 
mechanism  connected  directly  to  the  top  of  each  surface. 
Vanes  (1)  and  (3)  are  operated  together  as  " elevators n  and 
impart  a  moment  about  the  Y-axis  (pitch) .  Vanes  (2)  and  (4) 
together  are  the  " rudder n  and  contribute  a  moment  about  the 
Z-axis  (yaw).  Vanes  (1)  and  (3)  displaced  in  opposite 
directions  and  (2)  and  (4)  displaced  oppositely  work  as 
"ailerons"  to  impart  a  moment  about  the  X-axis  (roll) . 
Throughout  this  work,  the  term  aileron,  rudder,  and  elevator 
will  be  used  as  defined  here.  The  actual  torque  applied  by 
each  combination  of  vanes  was  determined  experimentally  and 
"constants  of  effectiveness"  were  calculated  by  the  AROD 
project  engineers. 
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The  constants  of  effectiveness  are  referred  to 
as  dimensional  stability  derivatives  when  developed  theoret¬ 
ically.  They  are  given  the  symbols  L,#,  Ntc,  Nr ,  ,  for  their 
contribution  of  moments  about  the  roll,  pitch,  and  yaw  axes 
due  to  a  displacement  by  the  ailerons,  elevators,  and 
rudder.  They  can  be  thought  of  as  acceleration  quantities 
per  radian  of  control  surface  displacement.  The  resulting 
relationships,  concluding  in  moments  about  the  three  body 
axes,  are 
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(2.21) 


**  -  I*.*. 

*t  -  H.,4* 

Ht  -  NP,Sr 

where  s  is  the  displacement  of  the  aileron ,  elevator 
or  rudder  into  the  rotor  downwash.  L,#,  M,#,  and  NP, 
are  given  in  Appendix  A. 

Theoretically,  each  set  of  control  surfaces  will'  yield  a 
constant  of  effectiveness  about  each  axis  (three  constants 
for  each  aileron,  elevator,  and  rudder) .  The  symmetry  of 
AROD  causes  negligible  cross  coupling  of  control  surfaces, 
however,  and  it  is  ignored.  Note  also  that  the  effective¬ 
ness  of  each  control  vane  was  determined  for  one  rotor  speed 
and  that  the  change  in  downwash  associated  with  the  changes 
in  rotor  speeds  is  also  ignored. 

3.  Servo  Equations  for  Control  Surfaces  and  Throttle 
Commanded  inputs  sent  to  AROD  include  "go  for¬ 
ward/backward",  "go  side  to  side",  "turn  around",  and  "go 
up/down".  These  commands  can  be  satisfied  by  pitching  the 
AROD  forward  or  sideways  (pitch  and  yaw)  for  translational 
flight,  turning  the  AROD  about  its  centerline  (roll)  for 
turning  about,  and  changing  the  propeller  speed  for  altitude 
control  (throttle) .  second  order  servo  motors  are  used  to 
displace  each  control  vane  as  well  as  the  throttle, 
a.  Control  Surface  Servos 

The  servos  for  each  control  surface  are  identi¬ 
cal.  Three  equations  will  describe  the  operation  of  the 
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elevators,  rudders  snd  ailerons.  For  each  of  thasa  aqua¬ 
tions,  at  least  two  servos  are  actually  operating  at  the 
same  tine  on  different  control  vanes.  Bach  servo  receives  a 
command  input  (voltage)  and  results  in  an  angular  position 
for  the  servo. 

4.  -  -Hi 4,  -  H2 5,  +  HjU.  (2.22) 

4a  ■  -Hi 4a  -  H24,  +  HjU, 
fir  “  *Hifir  -  H|4P  +  HjU,. 

where 

H1  and  H2  are  servo  gain  constants, 
u, ,u# ,  and  ur  are  the  servo  input. 

b.  Throttle  Servos 

The  servo  motor  used  to  open  and  close  the 
throttle  is  identical  to  the  ones  used  for  the  control 
surfaces.  An  additional  first  order  system  describes  the 
relationship  between  the  throttle  position  and  the  power 
delivered  to  the  rotor  (thrust,  TP ) .  The  complete  third 
order  system  is 

St  m  — Hi 4^  —  H24t  +  H2 Uf  (2.23) 

Tr  -  -Tr  +  .  5  (4t  ) 

The  power  delivered  to  the  rotor,  TP  ,  directly  relates  to 
both  the  change  in  rotor  speed,  AwP ,  and  the  force  in  the  X 
direction.  These  relationships  are 

A«r  -  TP  *  Pr#f  (2.24) 

AFt„-  TP 

where  Prcf  is  a  constant  relating  the  rotor  speed  to  the 
power  delivered  by  the  engine. 
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D.  ASSEMBLING  THE  MODEL:  PROGRAMMING 


In  this  section,  ths  squat iors  of  motion  and  the 
force/moment  relationships  are  represented  in  a  block 
diagram  of  the  AROD  model  developed  earlier,  and  the  flow 
chart  used  to  program  the  model  is  illustrated. 

1.  The  Complete  Block  Diagram 

Figure  2.7  is  a  block  diagram  of  the  complete  model. 
For  a  more  detailed  explanation  of  the  experimental  aero¬ 
dynamic  forces  and  moments ,  Appendices  A  and  B  should  be 
consulted. 

2.  Program  Flow  Chart 

The  logic  flow  chart  developed  to  program  the 
equations  and  the  aerodynamic  data  is  pictured  in  Figure 
2.8.  The  IBM  Dynamic  Simulation  Language  (DSL/VS)  source 
code  used  is  given  in  Appendix  B  [7]. 

E.  SIMULATION  RESULTS 

Graphical  illustration  of  time  response  is  shown  here. 
The  four  inputs  u, ,u* ,ur ,  and  u,  are  applied  alternately  as 
step  inputs  to  the  AROD  model .  A  caution  to  the  reader 
interpreting  these  results  is  appropriate  here.  The  body- 
fixed  angle  rates  P,  Q,  and  R  are  the  dynamic  states  of 
interest  in  the  control  problem.  Emphasis  should  be  placed 
on  the  response  of  these  rates  to  control  input  and  pot  on 
the  earth-fixed  Euler  angles  p,  6,  and  i>  .  The  Euler  angles 

are  shown  here  only  as  a  reference  for  the  reader  who  wishes 

• 
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Figure  2.8  AROD  Simulation  Logic  Flow  Chart 
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to  picture  what  the  vehicle's  behavior  would  be  relative  to 
the  Earth.  Confusion  can  result  if  one  tries  to  directly 
associate  these  angles  with  the  roll,  pitch,  and  yaw.  By 
way  of  example  (refer  to  Figure  2.1),  the  time  derivatives 
of  ip,  0,  and  are  equal  to  roll(P),  pitch(Q),  and  yaw(R) 
only  when  <p  -  0  and  0  and  1>  are  small.  If  the  body  has 
rotated  90°  ( <p  =  90),  then  Q  actually  approximates ~  the 
derivative  of  ^  and  R  the  derivative  of  0. 


Figure  2.9  is  the  resulting  time  response  of  the 
model  when  all  control  inputs  are  zero.  For  zero  input,  all 
control  surfaces  are  parallel  with  the  X-axis  and  not 
displaced  into  the  downwash.  The  throttle  is  set  for  hover 
speed.  Note  that  the  pitch  and  yaw  axes  remain  stable  and 
zero  as  expected,  but  the  roll  axis  spins  opposite  to  the 
rotation  of  the  propeller.  This  precession  is  a  result  of 
downwash  swirl  which  was  measured  and  included  in  the 
tabular  aerodynamic  data.  The  propeller  downwash  is  not  a 
linear  flow  of  air  but  has  a  pattern  of  vortices  which 
results  in  a  force  applied  to  the  control  surfaces  leading 
to  a  moment  about  the  roll  axis. 


Figure  2.10(a)  is  the  resulting  time  response  when 
the  aileron  (sa )  is  displaced  and  all  other  control  inputs 
are  zero  as  in  Figure  2.9.  Figure  2.10(a)  indicates  that 
the  ailerons  can  be  positioned  to  counteract  the  swirl 
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Figure  2.10  (a)  Time  Response  to  Aileron  Input 


effect  and  actually  roll  the  AROD  oppositely.  These  results 
also  affirm  that  the  roll  axis  Is  uncoupled  from  the  pitch 
and  yaw  axes  as  Equation  2.9  suggests.  Figure  2.10(b)  is  an 
illustration  intended  to  show  that  the  roll  can  be  con¬ 
trolled  and  set  to  zero  when  no  roll  is  desired.  In  this 
example,  an  initial  displacement  of  -.0499  on  the 
aileron  control  vanes  is  used  and  a  bias  equal  to  this 
initial  displacement  is  added  to  ua . 

3.  Elevator  and  Rudder  Input 

Figure  2.11  is  the  result  when  the  rudder  combina¬ 
tion  of  vanes  (Sr  )  is  displaced  into  the  downwash  and  the 
ailerons  are  set  as  in  Figure  2.10(b)  to  minimize  roll. 
Note  the  severe  coupling  between  the  pitch  and  yaw  axes,  as 
expected.  This  coupling  is  a  major  complication  to  the 
control  process  and  is  dealt  with  in  the  next  chapter. 
Figure  2.12  is  the  result  when  the  elevator  combination  of 
vanes  (S# )  is  displaced  into  the  downwash,  rudder  is  reset 
to  zero,  and  the  ailerons  set  as  above.  The  results  are 
similar  to  the  rudder  observations.  .  In  both  cases,  the  roll 
has  been  affected  by  the  pitch  and  yaw  disturbances,  a 
result  of  the  tabular  aerodynamic  data  which  changes  with  a 
and  fhe  amplitude  of  the  elevator,  however,  is  greater 
than  the  illustration  with  the  rudder  and  more  clearly  shows 
the  serious  problems  that  these  aerodynamic  moments  present. 
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Figure  2.10  (b)  Time  Response:  Aileron  Displaced 
to  Eliminate  Swirl  Effect 
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Figure  2.12  Time  Response  to  Elevator  Displacement 
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Figure  2.13  is  the  result  when  the  throttle  ( it )  is 
opened  to  increase  the  rotor  speed  and  all  other  inputs  are 
zero.  The  rotor  applies  a  torque  about  the  roll  axis  and, 
indeed,  this  influence  is  seen  in  the  results.  The  torque 
is  directed  oppositely  to  the  spin  direction  (-P)  and 
results  in  a  less  negative  roll.  However,  the  equations'  for 
pitch  (Q)  and  yaw  (R)  include  terms  coupled  to  the  propeller 
speed  through  the  aerodynamic  moments.  The  aerodynamic 
moments  are  determined  by  angle  of  attack  which  is  a 
relationship  between  the  body- fixed  velocities.  As  U 
increases,  the  coupling  with  pitch  and  yaw  becomes  more 
pronounced,  the  angle  of  attack  changes,  and  the  moments 
determined  from  the  tables  change  radically.  The  discon¬ 
tinuity  is  best  observed  by  the  sharp  change  in  roll,  pitch 
and  yaw  at  .75  seconds  in  Figure  2.13. 


It  is  insightful  to  observe  the  behavior  of  the 
model  with  the  troublesome  aerodynamic  moments  removed.  The 
system  behavior  without  the  wind  tunnel  data  supplying  the 
moments  is  much  more  in  keeping  with  expectations.  Figure 
2.14  is  the  response  of  this  modified  system  to  a  rudder 
displacement.  Note  the  similarities  of  the  nutation  (Q  and 
R)  and  precession  (i>)  to  the  gyroscope  illustration  of 
Figure  2.2.  Figure  2.15  is  the  response  of  the  same  system 
to  an  opening  of  engine  throttle.  This  smooth  system 
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Figure  2.14  Time  Response  to  Rudder  Control 
Aerodynamic  Moments  Removed 
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Figure  2.15  Tine  Response  to  Throttle  Input 
Aerodynamic  Moments  Removed 
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response  is  clearly  mors  in  kssping  with  how  ws  expect  ths 
actual  system  to  behave. 

The  validity  of  the  data  was  not  able  to  be  con¬ 
firmed  nor  is  it  being  challenged  here.  However,  these 
results  do  suggest  that  the  application  of  the  data  to  the 
system  model  is  a  problem  that  is  best  addressed  by  the 
engineers  who  derive  it. 
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The  equations  of  motion  describing  AROD  and  the  forces 
and  moments  which  act  on  AROD  were  developed  in  detail  in 
Chapter*  II.  A  computer  simulation  of  these  equations 
indicated  solid  stability  in  hover,  but  an  inability  to 
remain  stable  once  disturbed  from  the  hover  state.  There¬ 


fore,  a  method  is  needed  to  maintain 


and  drive  the 


system  to  a  steady-state  condition.  Additionally,  once  a 
steady-state  is  achieved,  it  is  necessary  to  steer  the  AROD 
up/down,  f orward/baclcward ,  side  to  side,  and  make  it  turn 
.about  its  axis  (roll)  .  Hence,  a  method  is  also  needed  to 
control  tie  AROD  and  cause  it  to  change  from  a  given  initial 
steady-state  to  another  desired  steady-state.  Stability  and 
control  are  the  topics  of  this  chapter.  Specifically,  an 
algorithm  based  on  optimal  control  theory  is  developed  in 
the  form  of  the  regulator  problem  to  address  the  stability 
issue.  The  regulator  is  then  adapted  to  the  control  issue 
formulated  as  the  tracking  problem  of  control  theory. 
First,  however,  a  summary  of  the  entire  stability  and 
control  issue  is  offered.  Then,  a  simplification  of  the 
complete  AROD  system  is  made  so  that  standard,  well  founded 
methods  may  be  used  to  address  the  stability  and  control 
dilemma. 
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A.  SUMMARY  OF  THE  CONTROL  PROBLEM 


The  task  is  to  provide  an  algorithm  on  which  to  base  a 


control  system  for  AROD.  The  control  system  intended  for 


AROD  relies  on  inputs  provided  by  an  operator  to  command 


desired  outputs.  This  section  discusses  the  control  problem 


based  on  the  inputs,  the  outputs,  and  the  equations  which 


relate  them. 


The  result  of  Chapter  II  was  18.  first  order  dif¬ 


ferential  equations  which  are  listed  in  Table  3.1.  Also 


listed  are  the  equations  for  useful  quantities  which  were 


derived  in  modelling  the  system. 


Figure  3.1  is  a  simple,  open  loop  block  diagram  of 


the  AROD  system.  The  block  labelled  "AROD"  represents  all 


Figure  3 . 1  AROD  Open  Loop  Model 
the  equations  listed  in  Table  3.1.  The  forces,  including 
gravity,  and  moments  acting  on  the  block  are  external 
influences  and  are  variables  in  the  equations  over  which  the 


designer  has  no  control.  The  inputs  of  interest  are  uB ,  ue , 


2.  Optimal  Feedback  Control 


Control  systems  may  be  divided  into  two  basic 
categories:  open  loop  and  closed  loop  systems.  Open  loop 

control  was  attempted  in  Chapter  II  and  is  pictured  in 
Figure  3.1.  Step  inputs  for  pitch (Q)  and  yaw(R)  resulted  in 

severe  nutation  and  oscillations.  Open  loop  control  is  not 

/ 

suited ‘  to  the  AROD  application.  Closed  loop  (feedback) 
control,  then,  is  the  direction  taken  in  this  work. 

A  feedback  control  system  is  one  which  tends  to 
maintain  a  controlled  variable  at  a  fixed  steady-state  value 
[8].  Therefore,  a  control  system  must  be  able  to  sense 
variations  in  the  system  variables  of  interest  in  order  to 
apply  control  to  them.  For  AROD,  the  outputs  are  the 
controlled  system  variables.  A  scheme  to  relate  the  output 
to  the  input  is  the  feedback  loop  and  is  pictured  in  Figure 
3.2.  The  feedback  matrix,  K,  must  be  determined  so  that  a 
given  control  input  always  results  in  the  same  desired 
output.  Note  that  the  reference  control  vector,  rc  ,  does 
not  match  the  commanded  input  vector,  uc  ,  except  in  steady- 
state.  The  time  when  rc#  u,.  is  called  the  transient 
response  time  and  is  also  the  time  it  takes  the  output  to 
reach  the  desired  state.  It  follows  that  transient  response 
time  should  be  as  short  as  possible.  Note  also  that  the 
difference,  rc -uc ,  should  be  greatest  when  rc  is  first 
applied  and  should  decay  to  zero  as  the  outputs  reach  their 
desired  state.  This  will  result  in  a  smooth  transition  of 
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Figure  3.2  Output  Feedback  Model 
the  output  states  from  beginning  to  end.  In  other  words, 
there  should  be  little  or  no  overshoot  of  the  desired  state 

by  the _ transient  states.  A  good  feedback  control  will 

minimize  both  the  transient  response  time  and  the  overshoot 
of  the  desired  state. 

A  method  of  feedback  control  that  might  minimize 
both  the  transient  response  time  and  the  overshoot  is 
optimal  control.  Its  name  suggests  that  optimal  control  may 
provide  the  optimum  solution  to  the  problem  of  feedback 
control  for  AROD.  However,  it  is  more  than  semantics  which 
drives  the  choice  of  controllers. 

First,  optimal  control  lends  itself  nicely  to  a 
discrete  time  solution  of  the  control  problem.  AROD  will 


employ  an  on  board,  digital  computer  to  perform  inflight 
stability  and  control.  While  a  continuous  time  controller 
can  be  easily  discretized  in  many  cases,  design  of  a  sampled 
data  controller  will  simplify  the  procedure. 

Second,  optimal  control  is  ideally  suited  to  state 
variable  techniques.  A  state-space  approach  is  desired 
because  of  the  enormity  of  the  system  and  the  complexity  of 
equations.  Matrix  manipulation  utilizing  computers  will 
prove  far  less  tedious  than  traditional  transfer  function 
methods . 

Third,  and  by  far  the  most  important  reason  to 
select  optimal  control  analysis  over  classical  methods  is 
the  unique  ability  of  optimal  control  theory  to  accommodate 
multiple-input,  multiple-output  systems  such  as  AROD  [3]. 
Four  inputs  and  four  outputs  have  been  identified.  It  will 
be  shown  in  the  next  section  that,  although  some  separation 
of  states  is  possible,  severe  coupling  between  other  states 
precludes  a  reduction  of  the  system  to  single-input,  single¬ 
output  systems . 

Fourth  and  last,  optimal  control  is  known  to  provide 
robust  and  insensitive  solutions  to  the  feedback  control 
problem.  Assuming  that  an  appropriate  performance  measure 
is  chosen  to  determine  the  optimal  feedback  gain  matrix,  K, 
the  solution  can  be  expected  to  have  a  fair  degree  of 
tolerance  to  plant  model  inaccuracies.  This  section 
develops  an  approximation  of  the  full  AROD  model  so  that 
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troublesome  nonlinearities  can  be  dealt  with.  Therefore , 
some  inaccuracies  are  bound  to  appear  in  the  new,  linear 
model  compared  with  the  nonlinear  system  of  Table  3.1. 
Clearly,  robustness  is  not  only  a  desired  property  of  the 
controller,  it  is  an  absolute  necessity  if  the  controller  is 
to  be  applicable  to  both  the  linear  and  nonlinear  model  of 
AROD.  *-  ~  ' 

With  the  advantages  of  optimal  control  comes  a 
certain  drawback:  optimal  control  requires  the  availability 
of  all  system  states  [3:p.  22] .  Therefore,  the  feedback 

scheme  pictured  in  Figure  3.2  must  be  modified  to  include 
all  system  states  in  the  feedback  loop.  Figure  3.3  satis¬ 
fies  this  requirement.  However,  if  the  system  states 


Figure  3.3  Full  State  Feedback  Model 
are  not  available  from  sensor  measurements  or  other  means, 
computational  observers  must  also  be  included  to  provide  an 


estimate  of  the  hidden  states.  Observers  add  an  extra 
computational  task  to  the  onboard  computer,  but  otherwise 
present  no  problem  to  the  control  designer. 


B.  LINEARIZING  THE  AROD  MODEL 

Throughout  model  development  in  Chapter  II,  nonlinear 
relationships  were  identified.  These  nonlinearities  fall 
into  two  categories:  nonlinear  combination  of  states  (e.g. 
Equation  2.8)  and  discontinuous  functions  (e.g.  table  look¬ 
up  of  aerodynamic  force/moment  coefficients) . 

A  linear  approach  to  controller  design  is  much  preferred 
over  a  nonlinear  analysis.  Three  major  reasons  drive  this 
preference.  First,  the  theoretical  basis  for  controllers  of 
linear  systems  is  well  developed  and  easier  to  implement. 
Second,  in  airplane  design,  the  use  of  linear  approximations 
for  nonlinear  systems  has  yielded  good  results  for  the  case 
of  steady-state  flight  and  small  disturbances  theref rom  [ 5 ] . 
Third,  the  optimal  control  solution,  K,  can  often  be  reduced 
to  a  single,  constant  gain  schedule  for  linear,  time- 
invariant  systems.  A  constant  gain  matrix,  K,  will  mean 
that  a  minimum  amount  of  memory  storage  will  be  required  of 
the  on  board  computer  for  control  implementation  [9]. 

In  this  section,  the  equations  of  Table  3.1  will  be 
replaced  with  linear  approximations  based  on  the  small- 
disturbance  theory,  steady-state  assumptions,  and  selected 
physical  approximations  where  necessary. 


1.  The  Small-Disturbance  Theory 


A  straight  forward  discussion  of  the  small  distur¬ 
bance  theory  is  contained  in  [5:p.  106].  Essentially,  if 
AROD  is  assumed  to  be  in  a  certain  steady-state  reference 
condition,  then  motion  of  the  vehicle  consists  of  small 
deviations  from  the  reference.  If  this  steady-state  is 
hover  fall  translational  and  angular  movement  is  aero)  “then 
each  of  the  states  in  Table  3.1  will  be  zero  or  nearly  zero. 
This  is  not  a  bad  assumption  in  the  case  of  AROD.  The 
flight  condition  in  which  AROD  will  spend  most  of  its  time 
is  the  hover.  Only  small  deviations  from  the  hover  will  be 
necessary  to  effect  translational  or  angular  movement.  For 
example,  a  translational  velocity  forward  (along  the  Z' 
axis)  of  four  miles  per  hour  requires  only  a  15  degree  pitch 
of  the  body.  The  small -disturbance  theory  and  the  steady- 
state  hove**  assumption  result  in  the  following  simplifica¬ 
tions  : 

(1)  Products  of  states  are  zero  since  the  product  of  two 
small  values  is  an  extremely  small  value. 

(2)  Derivatives  of  states  are  zero  since  the  time 
derivative  of  a  sma2 1  valued  time  function  is 
extremely  small. 

(3)  The  sine  of  a  state  is  equal  to  the  state  and  the 
cosine  of  a  state  is  equal  to  one.  This  is  the 
small  angle  approximation  for  angles  less  than  15 
degrees . 

Table  3.2  is  the  set  of  equations  which  result  when  the 
above  assumptions  are  applied  to  Table  3.1.  Note  that  much 
of  the  coupling  between  states  and  all  of  the  nonlinear 
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State 

• 

■  p 

• 

e 

-  Q 

*  .  • 

* 

-  R 

• 

p 

-  -(IP/IX  )wP  +  (L,  +  L*)/^ 

• 

Q 

-  -(Ir/Iy)wPR  +  (M,  +  Mt)/Iy 

• 

R 

-  (IP/IX)«PQ  +  (N.  +  Nt  )/I2 

• 

U 

-  g  +  Ftx/m  +  Ftx/m 

• 

V 

-  F.y/m  +  Fty/m 

w 

-  Fax/m  +  Ftz/m 

sa 

*  “H1 "  H25«  +  H2U. 

•• 

*• 

»  -H^*  -  Hz  5.  +  H2u, 

*; 

-  -  H2«P  +  H2ur 

•• 

• 

-  -Hj  fit  -  H2  5t  +  H2  ut 

products  of  states  have  been  eliminated.  Note  that  the 
derivative  term  of  rotor  velocity,  wP ,  remains  in  the 
equation  for  roll  (P) .  The  term  is  kept  because  u>r  is  not  a 
state  and  because  it  is  easily  approximated.  We  wish  to 
keep  as  much  similarity  between  the  linear  and  nonlinear 
models  as  possible,  so  physical  approximations  of 
troublesome  quantities  are  important. 
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2.  Phvalcal  Approximations  of  Forces  and  Moments 


The  force  and  moment  terms  in  Table  3.2  are  made  up, 
in  part,  of  functions  of  aerodynamic  coefficients  measured 
at  discrete  AROD  flight  attitudes  in  a  wind  tunnel.  The 
data,  listed  in  Appendix  A,  follow  no  linear  pattern  and, 
instead,  are  discontinuous  functions  of  angle  of  attack. 
The  relationships  used  to  include  these  data  in  '  the  non¬ 
linear  model  were  discussed  in  Chapter  II.  If  a  steady- 
state  condition  is  assumed,  then  some  simplification  of 
these  forces  and  moments  is  possible. 

It  is  not  desirable  to  exclude  the  aerodynamic  data 
completely  from  the  simplified  model.  An  analysis  of  the 
magnitude  of  the  forces  and  moments  which  result  reveals 
that  they  are  significant  in  all  flight  attitudes.  However, 
if  a  constant  hover  is  assumed  then  only  one  set  of  coeffi¬ 
cients  need  be  considered.  Hover  condition  is  defined  by 
AROD  project  engineers  as  a-90°  and  0-0°  [6],  In  hover 
condition,  the  forces  and  moments  become  those  listed  in 
Table  3.3.  Note  that,  for  the  purpose  of  Table  3.3,  the 
total  body-fixed  velocity  is  reduced  to  U.  V  and  w  are 
assumed  zero. 

Table  3 . 3  indicates  that  each  force  and  moment  in 
hover  flight  can  be  reduced  to  a  linear  function  of  a  state, 
with  one  exception.  Recall  from  Chapter  II,  we  concluded 
that  the  rotor  downwash  created  a  swirl  resulting  in  a  force 
on  the  control  vanes  and  a  moment  about  the  roll 
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TABLE  3.3 

1 

Forces 

Momenta 

AarQdttnamig 

Aerodynamic 

o 

L,-  U ( . 036) 

-  TP  (.045)  -  1.62 

F,y«  U (  . 0729) 

M,-  U( . 81) 

F,t«  -U(l. 08) 

N,-  -U  ( .  18) 

Thrust 

Thrust 

Ftx»  (76.5)  +  TP  (2.125)  L*  ■  -Tr(.5)  -  4,(21.29) 

o 

1 

>. 

Mt»  -4,  (14.51) 

Ft,»  0 

Useful  Approximations 
wr-  753.98  rad/sec 

Nt-  -4P  (16.68) 

«r-  Tr  (10.472) 

I 


axis.  This  moment  was  illustrated  in  Figure  2.10(a).  This 
same  constant  moment  appears  in  the  equation  for  roll 
moment,  I*  ,  in  Table  3.3.  The  compensating  factor  of  .0499, 
determined  experimentally  in  the  last  chapter,  can  be 
derived  by  dividing  the  constant  moment,  -1.06,  by  the 
aileron  control  vane  effectiveness,  L,,  (-21.29).  This 

constant  moment  will  be  eliminated  in  the  linear  model,  but 
wiJ  1  be  included  in  the  resulting  controller  as  a  constant 
addition  to  the  aileron  control,  u,  .  The  aerodynamic 
moments  about  the  pitch  and  yaw  axes  are  zero  in  hover 
flight.  While  this  poses  no  problem  to  the  linearizing 
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goal,  it  means  that  momants  in  othar  flight  attitudas  have 
not  baan  approximatad .  Thus  tha  controller  resulting  from 
the  linear  modal  will  not  consider  these  external  momants. 


Tha  complete,  linearized  equations  are  listed  in 
Table  3.4.  Numerical  values  for  the  constant  quantities 


were  taken  from  Appendix  A.  The  next  section  will  assemble 
the  equations  in  a  form  suitable  to  determine  the  optimal 
gain  matrix,  K. 


C.  FORMULATING  THE  CONTROL  PROBLEM 
1.  Stata  Space  Representation 

Tha  control  system  of  Figure  3.3  is  callad  a 
regulator.  As  discussed  earlier,  tha  goal  of  this  regulator 
aystam  is  to  obtain  a  stable  closed  loop  system  which  will 
force  any  arbitrary  initial  states  to  zero.  Tha  first  task, 
then,  Is  to  assemble  tha  stata  aquations  in  a  form  suitable 
for  control  analysis.  The  form  desired  is  given  by 

x(t)«  Ax(t)  +  BUe  (t)  (3.1) 

where 

x(t)  is  the  state  vector, 

A  is  the  plant  matrix, 

B  is  the  control  matrix, 
ue (t)  is  the  control  vector. 

Note  that  the  A  and  B  matrices  do  not  vary  with  time.  A 
time-invariant  system  greatly  simplifies  controller  design. 

The  18  states  listed  in  Table  3.4  are  the  complete 
states  of  the  equations  of  motion  for  AROD.  However,  some 
of  the  states  are  unnecessary  to  the  control  issue  and  some 
additional  assumptions  will  aid  in  determining  the  best 
control  scheme . 

The  control  problem  can  be  simplified  if  it  is 
recognized  that  translational  velocities  over  the  ground  (V' 
and  W'  in  the  earth-fixed  horizontal  plane)  are  accomplished 
by  pitching  the  AROD  forward/backward  or  side  to  side. 
Therefore,  the  states  important  to  this  task  are  roll(P), 
pitch (Q) ,  and  yaw(R)  along  with  their  associated  angles. 
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Tha  only  translational  valocity  of  intarast  is  tha  earth- 
fixad  vartical  valocity  or  altitude  rata.  For  small  pitch 
and  yaw  angles,  than,  body-fixad  V  and  W  are  equivalant  to 
earth-fixed  V*  and  W  and  will  ba  eliminated  from  tha 
aquations.  Also,  tha  earth-fixed  vartical  valocity  (al¬ 
titude  rata)  will  ba  approximated  by  tha  body-fixad  valoci¬ 
ty,  U.*  ~ 

A  caution  about  tha  anglas  ?,8,^,  is  appropriate 
here.  The  linearized  equations  conveniently  eliminated  any 
dependence  on  the  earth-fixed  Euler  angles.  While  this  may 
be  a  valid  assumption  for  determining  the  control  solution, 
it  changes  the  meaning  of  the  angle  states  6  and  <p 

depends  only  on  the  angles  6  and  which  will  be  kept  small 
in  the  hover  condition,  but  ©  and  i>  both  depend  on  the 
roll  (heading)  angle,  <p%  which  can  range  from  0  to  360°. 
Therefore,  the  linearized  angles  %  and  will  be  used  from 
here  on  to  remind  U3  that  these  are  body-fixed  pitch  and  yaw 
angles  and  no  longer  relate  the  vehicle  to  the  earth-fixed 
coordinate  system. 

The  complete,  linearized  system  can  now  be  put  in 
the  form  of  Eq.  3.1  and  is  listed  in  Table  3.5.  Figure  3.4 
is  the  complete  signal  flow  graph  for  the  linear  system. 

2.  Stability.  Observability,  and  Controllability 
a.  Stability 

A  gyroscope  is  a  stable,  open  loop  device. 
Small,  bounded  disturbances  of  the  roll,  pitch,  or  yaw  axes 


TABLE  3.5 


STATE-SPACE  UNEARIZED  SYSTEM 
X  (t) a  jjp  /6|j  jij)  b  ,U,P, Tr  fS^  6p  fS ^  r  &  $  t^r  /  3  ^ 

A(t)  = 


0 

0 

0 

0 

1 

0 

0 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

0- 

0 

0 

0 

0 

1 

0  0 

0  0 

0 

0 

t) 

o' 

0 

0 

0 

0 

0 

.894 

0 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

0 

0 

.024 

0  -3.74 

0 

0 

-21.29  0 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0  0 

0  .5 

0 

0 

0 

0 

0 

0 

0 

.483 

0 

0 

0  -6.75  0  -14.51 

0  0 

0 

0 

0 

0 

0 

0 

0 

.108 

0 

0 

6.78 

0 

0  0 

-16.68  0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0  0 

0  0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 

0  0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 

0  0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 

0  0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

-157.9  0 

0  0 

-17.77  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  -157. 

9 

0  0 

0 

-17.77  0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 

-157.9  0 

0 

0  - 

-17.77 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 

0  -157. 

9  0 

0 

0  -17 

.77 

“  0 

0  0 

0  0 

0 

0  0  0  0  0 

0 

157.9  0 

0 

0  " 

t 

B(t) 

= 

0 

0  0 

0  0 

0 

0  0  0  0  0 

0 

0  157. 

9  0 

0 

0 

0  0 

0  0 

0 

0  0  0  0  0 

0 

0  0 

157. 

9  0 

0 

0  0 

0  0 

0 

0  0  0  0  0 

0 

0  0 

0 

157.9 

Note:  t  is  the  transpose  operator. 


produce  bounded  outputs.  The  outputs  m ay  be  in  the  form  of 
nutation  or  oscillations,  but  it  will  not  grow  with  time 
once  the  disturbance  is  removed.  This  stability  was  illus¬ 
trated  for  AROD  in  Figure  2.14  where  the  yaw  axis  was 
disturbed.  The  results  compared  favorably  to  the  response 
of  a  gyroscope. 

The  open  loop  linearized  system  may  be  tested 
for  stability  also.  Since  it  is  linear  time- invariant,  the 


simplified  system  of  Table  3.5  is  stable  if  and  only  if  the 
eigenvalues  of  the  continuous  time  plant  matrix,  A,  have 
real  parts  all  less  than  zero  [10].  Our  system  fails  this 
test  for  stability  since  it  results  in  seven  out  of  16 
eigenvalues  equal  to  zero.  Therefore,  the  linear  model  has 
lost  some  of  the  characteristics  of  the  nonlinear  model. 
However,  the  control  algorithm  that  stabilizes  the  linear 
system  and  minimizes  the  time  response  and  overshoot  may 
still  be  applicable  to  the  AROD. 
b.  Observed!  ility 

Another  relationship  to  test  the  applicability 
of  the  linearized  plant  to  AROD  is  observability.  An 
observable  system  is  one  for  which  all  states  may  be 
determined  from  observations  of  the  output.  This  is  an 
important  concept  since  an  optimal  controller  will  be 
useless  if  any  system  states  are  not  known.  To  determine 
the  observability  of  our  linearized  plant,  the  outputs  which 
are  available  from  direct  measurement  must  be  known.  Refer¬ 
ence  2  indicates  that  three  single  axis  rate  gyros,  one 
vertical  rate  gyro,  a  magnetometer  (magnetic  heading)  and  a 
barometric  altimeter  will  be  available  on  AROD.  Control 
vane  position  based  on  a  history  of  the  input  is  also 
assumed  to  be  known.  Table  3.6  is  a  listing  of  the  states 
which  will  be  measured  and  the  device  used  for  each. 

The  output  of  the  linear,  time-invariant  system 

is  given  by 
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y(t)  -  Hx(t) 


(3.2) 


where 

y(t)  is  the  output  vector, 

H  is  the  measurement  matrix 

dimensioned  o  x  n, 
n  is  the  number  of  states, 

o  is  the  number  of  measurements. 


TABLE  3.6 

Magnetometer 

<p  (heading) 

Rate  Gyro 

U,  P,  Q,  R 

control  History  (Servo) 

5 a  i 

and  from  Table  3.6 


H  = 


1  0  0  0  0  0  0 
0  0  0  1  0  0  0 
0  0  0  0  1  0  0 
0  0  0  0  0  0  1 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 


0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
1  0  0  0  0 
0  10  0  0 
0  0  10  0 
0  0  0  1  0 
0  0  0  0  1 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


Observability  of  a  system  due  to  its  measured 

I 

A 

states  can  be  determined  from  the  observability  matrix,  W0 , 


given  by  [11] 

W0=  [H  „  HA  „  HA2  . HA"'1]  (3.3) 

where 

W0  is  a  n  x  mn  matrix, 
m  is  the  number  of  states  measured. 


The  rank  of  W0  (the  number  of  independent  columns)  must 
equal  n  if  all  of  the  states  are  observable  from  output 
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measurements.  For  the  16  state,  linearized  AROD  system  of 
Table  3.5,  W0  has  rank  10.  Although  a  control  solution 
might  still  be  found  assuming  all  states  were  available,  it 
would  not  be  a  realizable  controller.  Therefore,  a  dif¬ 
ferent  approach  must  be  taken. 

Figure  3.4  reveals  that  the  roll(P)  and  the 
throttle(Tr)  coupled  together  and  the  pitch (Q)  and  yaw(R) 
states  are  likewise  coupled  as  expected  from  the  gyroscope 
analogy.  However,  these  two  subsystems  are  coupled  only 
through  a  force/moment  relationship  dependent  on  the  total 
velocity  of  the  body,  Vtot.  In  the  last  section,  Vtot  was 
approximated  by  U.  A  further  simplification  assumes  that 
the  moment  on  the  pitch  and  yaw  axis  due  to  vertical  speed, 
U,  is  negligible.  The  complete  16  state  linear  system  can 
now  be  considered  two  independent  subsystems  of  eight  states 
each  as  pictured  in  Figure  3.5.  State  space  representation 
of  the  two  subsystems  is  listed  in  Table  3.7. 

Each  subsystem  has  an  observability  matrix  of 
rank  eight.  Therefore,  the  states  required  in  the  subsystem 
configuration  can  be  obtained  from  the  measurements, 
c.  Controllability 

A  final  test  of  the  linearized  model,  before  any 
effort  is  expended  in  actual  controller  design,  is  to  deter¬ 
mine  whether  it  is  even  possible  to  find  a  control  sequence 
such  that  a  desired  final  state  may  be  reached.  This  con- 
\s  known  as  controllability  [lOsp.  205].  The  control- 


lability  matrix  is  given  by 


where 


Wc=  [B,  AB,  AZB,  ...  ,  A" - 1 B] 


Wc  is  a  nxin  matrix, 
i  is  the  number  of  control  inputs. 


(3.4) 


The  measure  of  controllability  is  the  rank  of  Wc  .  For  each 
subsystem,  the  rank  is  eight.  Hence,  a  control  sequence  can 
be  found  for  the  linearized  system. 


3.  Discretizing  the  System 
a.  Sampling  Frequency 

We  wish  to  obtain  a  discrete  time,  digital 
control  algorithm.  Therefore,  the  system  of  Table  3.7  must 
be  converted  to  a  sampled  system.  A  basis  for  determining 
the  sampling  interval,  AT,  is  needed.  The  sampling  theorem, 
based  on  the  Nyquist  criterion  requires  that  [11]:  '  ~ 

f,=  ^  >  2f*  (3.5) 

where 

fs  is  the  sampling  frequency, 

AT  is  the  sampling  period, 

fM  is  the  natural  frequency  of  the  system. 

Multiple-input,  multiple-output  systems  do  not 
have  a  single  natural  frequency.  However,  if  the  system  is 
considered  a  set  of  several  single-input,  single-output 
systems  then  a  number  of  different  bandwidths  are  observed. 

T1  e  highest  bandwidth  can  be  used  as  the  natural  frequency 

of  the  entire  system  for  sampling  purposes.  In  this  way, 
four  independent  systems,  with  coupled  terms  set  to  zero, 
are  extracted.  The  frequency  domain  transfer  functions  are 
listed  in  Table  3.8  along  with  their  bandwidths.  A  Bode 
analysis  yields  the  frequencies  listed.  Each  of  the 
transfer  functions  also  has  the  describing  equation  of  the 
second  order  servo  in  the  denominator.  Second  order 
analysis  yields  a  servo  natural  frequency  of  2  Hz.  since 
the  servo  frequency  is  the  highest,  it  will  be  considered 


TABLE  3.8 

Natural 

Frequency 

Roll 

<&(  s) 

La  •  H2 

.74  hz 

U,  (s)“ 

S2  (S2  +H,  S+H2  ) 

Fitch 

8(S) 

M*  •  H2 

.61  hz 

u.  (s)“ 

s2  ( S2  +H,  s+H2  ) 

Yaw 

\b(  S) 

Nr  «H2 

.65  hz 

Ur  (S) 

S2  (S2  +H,  S+H2  ) 

u 

Ufs) 

-.89H2 

.03  hz 

ut  (s) 

s(s2+H,  s+H2  )  (s+1) 

the  natural  frequency,  fN ,  of  the  multiple  system.  The 
sampling  frequency  is  now  given  as 

f,  >  2f„  «  4  HZ  (3.6) 

Equation  3. 6.  is  a  lower  bound  on  the  sampling 
frequency  that  will  allow  the  continuous  time  system  to  be 
reconstructed  once  sampled.  Selection  of  the  actual 
sampling  rate  is  based  on  engineering  judgment  which  calls 
for  a  factor  of  10  between  the  highest  frequency  component 
and  the  sampling  frequency  [12].  A  rate  of  25  Hz  meets  this 
criteria  and  yields  a  sampling  period  of  AT*. 04  seconds, 
b.  Analog  to  Digital  Conversion 

The  process  of  conversion  is  illustrated  in 
Figure  3.6.  The  block  labelled  D-A  is  a  digital  to  analog 
signal  converter.  The  block  labelled  A-D  is  an  analog  to 
digital  signal  converter.  The  AROD  block  is  the  nonlinear 

76 


i  ut-i  u*  u *  un  uw «  tat  urn  uwA/tfi/w  \3wuwswvas\ai 


Figure  3 . 6  Digital  System  feedback  Scheme 

model  which  has  been  linearized  here.  The  next  step  is  to 
combine  the  D-A,  linear  AROD  model,  and  the  A-D  to  obtain  a 
discrete  time  representation  of  the  linear,  time- invariant 
system.  since  the  present  linear  model,  consists  of  two 
subsystems,  the  task  is  to  convert  the  two  systems  and 
determine  two  control  systems  separately.  The  discrete 
version  of  the  closed  loop  system  is  pictured  in  Figure  3.7. 
The  discrete  form  of  the  system  Equation  3.1  is  given  by 

x(k+l)-  <f>x(k)  +  TUg  (k)  (3.7) 

k  ■  kAT, 
k+1-  kAT+AT, 

AT  is  the  sampling  period, 
q)  is  discretized  A  plant  matrix, 
r  is  discretized  B  control  distribution. 


where 


H  p{>  yoo 


Figure  3.7  Discretized  AROD  Feedback 


The  discrete  plant  and  control  distribution  matrices  are 
given  by  [11] 

$  m  eAAT  (3.8) 

*  ~  r  -  dsB 

where 

e  is  the  natural  logarithm  operator, 

s  is  the  Laplace  operator, 

ds  is  the  derivative  with  respect  to  s. 

Equation  3.8  can  be  solved  by  using  a  series 
expansion  of  the  matrix  exponential  and  integrated  using  a 
numerical  scheme  on  a  digital  computer.  The  program  called 
ORCONV  [13]  is  used  to  obtain  the  discrete  time  system 
matrices  of  Table  3.9. 


TABUS  3.9 


•  i  ki* 


Lgfflasrag 


4>PT  “  rPT  “ 


1  o  .0400  -.0030  -.0167  0  -.0002  0 

-.0003  0 

0  1  0  .  0350  0  .  0004  0  0 

0  .00001 

0  .0009  1  -.1466  -.8217  -.0015  -.0134  0 

-.0299  -.00003 

00  0  .9608  0  .0189  0  .0003 

0  .0007, 

0  0  0  0  .9010  0  .0275  0 

.0990  '  0 

0  0  0  0  0  .  9010  0  .  0275 

0  .0990 

0  0  0  0  -4.334  0  .4132  0 

4.3.4  0 

0  0  0  0  0  -4.334  0  .4132 

0  4.334  _ 

<t>QR  " 

ro*  * 

1  0  .0395  -.0054  -.0113  .0012  -.0001  .00001 

-.0002  .00001 

0  1  .0054  .0395  -.0010  -.0130  -.00001  -.0002 

-.00001  -.0002 

0  0  .9636  -.2667  -.5530  .0879  -.0090  .0010 

-.0203  .0016 

0  0  .2681  .9636  -.0768  -.6357  -.0009  -.0104 

-.0014  -.0234 

0  0  0  0  .9010  0  .  0275  0 

.0990  0 

0  0  0  0  0  .9010  0  .0275 

0  .0990 

0  0  0  0  -4.334  0  .4132  0 

• 

4.334  0 

0  0  C  0  0  -4.334  0  .4132 

»  m 

0  4.334 

D.  THE  OPTIMAL  REGULATOR 

A  review  of  the  linear  control  problem  is  appropriate 
here.  We  now  have  two  discrete-time  subsystems,  each 
defined  by  Eq.  3.7.  The  roll/throttle  subsystem  is  descri¬ 
bed  by  the  plant  matrix  <j>p  T  and  control  distribution  rp  T  . 
The  pitch/yaw  subsystem  is  described  by  4>QR  and  rQR  These 
matrices  were  listed  in  Table  3.9.  We  also  know  that  all  of 
the  states  in  each  subsystem  are  available  either  by  direct 
measurement  or  computational  observers.  We  wish  to  find  the 
closed  loop  feedback  matrix,  K.,  in  Figure  3.7  such  that  any 
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disturbance  of  any  atata  from  the  reference  hover  condition 
(all  states  zero)  will  be  damped  out,  returning  the  states 
to  zero  in  minimum  time  with  minimum  overshoot.  Since  each 
subsystem  is  controllable,  the  control  law  that  satisfies 
these  requirements  is  given  by  [10:p.  771] 

U«*(k)  -  -K.X(k)  .  (3*9) 

where  *  denotes  optimal  in  the  minimum  time,  minimum 
overshoot  sense. 

Note  that  K  is  a  constant  matrix.  We  state  this  goal  to 
reduce  to  a  minimum  the  amount  of  computer  memory  required 
to  store  the  feedback  gains  and,  in  general,  to  simplify  the 
control  algorithm.  The  validity  of  this  approach  is  born 
out  by  experience  with  linear  time- invariant  systems  [3]. 
The  applicability  of  the  resulting  controller  to  the 
nonlinear  system,  like  the  previous  linearizing  assumptions, 
will  be  determined  by  the  time  response  of  the  full  AROD 
model  once  K  is  found.  The  following  paragraphs  will 
provide  analytical  justification  for  the  steady-state 
approach . 

l-  Qptimai  gontyoi 

Optimal  control  theory  is  concerned  with  minimizing 
a  performance  criterion  chosen  by  the  designer  [3:p.  3]. 

Therefore,  the  primary  focus  in  developing  the  controller 
should  be  on  a  performance  measure  that  quantifies  the  goals 
of  the  design  effort.  These  goals  are  summarized  here: 

(1)  Minimize  the  transient  response  time 
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(2)  Minimize  the  state  overshoot 

(3)  Determine  a  constant  gain  schedule,  ft. 

(4)  Operate  within  the  physical  constraints  of  the 
system. 

The  performance  measure  typically  used  in  formulating  the 
discrete  time  optimal  regulator  control  solution  is 


J-  £x(H)  *&(]!)  (3.10) 

+  2*  E  [x(k)*Q(k)x(k)  +  u,.  W'RClOUe  (k)  ] 

k  »  0 

where 

J  is  the  cost  function, 

N  is  the  final  time  step  (NAT) , 
k  is  the  time  step  index, 

Q(k)  is  the  state  weighting  matrix, 

R(fc)  is  the  control  cost  weighting  matrix, 

C  is  the  final  state  weighting  matrix, 
t  is  the  matrix  transpose  operator. 


The  proper  choice  of  C,  Q(k),  and  R(k)  will  ensure  that  the 
control  schedule,  uc (k) ,  that  minimizes  J  will  meet  the 
first  two  of  our  design  goals  and  minimize  the  control 
effort  required. 

The  solution  to  Eq.  3.10  for  time-inva?*iant  systems 
is  found  from  the  recurrence  relation 

ft  (N-k)  -  -[R(N-k)  +  r^k-ljr]-1  (3.11) 

x  r*p(k-i)(() 

and 

P(k)*  M  (N-k)  *  P  (k-1)  M  (N-k)  +  K  (N-k)  *  R  (N-k)  K  (N-k)  +  Q  (N-k) 

and 

M(k)  =  (j)  -  nc  (N-k) 


To  meet  goal  ( 3 ) ,  a  constant  K,  It  is  necessary  to 
know  that  under  certain  conditions  the  recurrence  relation 
converges  to  a  steady-state  for  large  enough  k.  The 
conditions  sufficient  for  a  steady-state  convergence  are: 
constant,  positive  definite  weighting  matrices  Q  and  R; 
final  state  weight  OO?  and  uc (k)  unconstrained  [10].  The 
final  time  loses  its  meaning  for  the  steady-state  case  since 
we  assume  convergence  prior  to  reaching  time-step  N. 
Equation  3 . 10  becomes 

k  c 

J-  £  [x(k)tQx(k>  +  uc  (k) 1  Ruc  (k)  ]  (3.12) 

k  =  0 

where  kc  is  the  time  step  when  J  converges. 

The  recurrence  formula  can  now  be  written  [13] 


and 

and 


k  -  [R  +  r'pr]'  Vpt))  (3.13) 

P  =  M*  PM  +  K*  RK  +  Hl  QH 
M  =  (j)  -  TK 


where  the  measurement  matrix,  H,  modifies  Q  so  that  only 
states  of  interest  are  considered. 


The  revised  performance  measure  appears  to  have  no 
dependence  on  the  final  state.  In  fact,  the  solution,  K,  is 
sometime  referred  to  as  a  " subopt imal"  gain  schedule. 
However,  experience  shows  that  for  time-invariant  systems, 
the  suboptimal  gain  schedule  is  often  quite  satisfactory. [9] 
2.  Selecting  the  Performance  Measure 

The  weighting  matrices  Q  and  R  must  be  selected  to 
define  the  cost  function  used  to  determine  the  optimal 
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gains.  For  the  linear  AROD  problem,  Q  and  R  are  both 
diagonal  matrices.  Each  diagonal  element  of  Q  corresponds 
to  a  weight  placed  on  the  deviation  of  a  particular  state 
from  its  desired  value.  Equation  3.13  introduced  the  idea 
of  weighting  only  certain  states  utilizing  the  H  matrix. 
Tvje  states  we  are  interested  in  weighting  are  the  angle 
states,  the  altitude  rate,  and  the  control  vane  displace¬ 
ments  and  displacement  rates.  The  angle  states  and  altitude 
rates  describe  the  steady-state  of  the  system.  The  control 
vanes  are  constrained  in  their  movement.  Therefore,  these 
are  the  states  with  which  we  are  concerned  with  controlling 
through  the  gain  schedule.  A  definition  of  the  measurement 
matrix,  H,  strictly  for  the  purpose  of  the  cost  function,  is 

1000000  0]  (3.14) 

01000000 
H  =  00001000 

0  0000100 
00000010 

-  -  0  000000  1 
where  H  applies  to  both  subsystems.  The  non-zero 
elements  correspond  to  <p,TJ,sa  ,St  ,Sa  ,5*t  for  roll/ throttle 
and  0b  ,r/}b  ,$P  #*•  /S*r  for  pitch/yaw. 

Now  Q  is  a  6x6  matrix.  If  all  the  states  of  interest  are 
weighted  equally,  Q  is  given  by 

1  0  0  0  0  0]  (3.15) 

0  1  0  0  0  0 
Q  =  001000 

0  0  0  1  0  0 
0  0  0  0  1  0 
0  0  0  0  0  1 
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This  reduction  in  cost  states  will  streamline  the  iterative 
design  process  of  the  gain  matrix  in  the  next  section. 

Each  subsystem  has  two  control  inputs:  u,  and  ut 
for  roll/throttle,  and  u# ,  uP  for  pitch/yaw.  A  control 
weight  matrix  which  treats  each  control  cost  equally  is 

„  „  ft  °1 
I®  ij 

The  cost  function  resulting  in  the  best  gain  matrix 
for  our  purposes  may  not  be  the  one  which  weights  each  state 
and  cost  equally.  Factors  such  as  the  units  a  particular 
state  is  scaled  to  and  system  constraints  on  other  states 
require  scaling  the  diagonal  values  of  the  Q  and  R  matrices 
accordingly.  It  may  also  be  more  important  to  drive  a 
particular  state  more  rigidly  than  others  .  because  of  the 
dynamics  of  the  system.  In  the  case  of  AROD,  it  is  most 
important  to  drive  the  angle  states  and  the  altitude  rate  to 
the  desired  state  because  the  other  states,  all  derivatives 
of  these,  follow.  The  time  response  of  the  angle  states  and 
altitude  rate  is  the  standard  by  which  the  best  cost 
function  is  selected.  Therefore,  an  iterative  approach  is 
taken,  altering  Q  and  R  as  necessary  to  meet  the  following 
specific  objectives:  settling  time  of  two  seconds  to  less 
than  10%  overshoot,  and  maintaining  constrained  states 
within  limitations. 

The  settling  time  requirement  is  chosen  for  two 
reasons.  First,  the  servo  response  time  (second  order  time 
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constant)  is  .5  seconds  which  places  a  lower  limit  on  the 
response.  The  equations  are  coupled  which  means  action  by 
one  state  will  force  a  delayed,  but  proportional  reaction  by 
the  others.  Two  seconds  for  any  response  on  a  slow  speed 
hovering  vehicle  such  as  AROD  is  considered  adequate. 
Second,  if  the  erratic  behavior  observed  in  Chapter  II  is  to 
be  controlled,  then  a  quick  response  to  disturbances  and 
commands  is  necessary. 

The  control  vanes  are  constrained  in  displacement  to 
+.5236  radians  (30°)  and  in  velocity  to  +.8727  radians  per 
second.  The  throttle  is  constrained  to  +100  radians  and 
+100  radians  per  second.  The  iterative  design  process  must 
account  for  these  constraints  by  weighting  the  corresponding 
elements  of  the  Q  matrix.  Appendix  A  lists  constraints. 

Constraints  on  the  servo  states  imply  that  the  rate 
of  change  of  other  states  may  be  restricted  as  well.  A 
constraint  on  the  output  states  means  a  limit  must  be 
applied  to  the  disturbances  from  which  we  expect  the  system 
to.  recover  in  the  given  time.  An  angular  displacement  of  .1 
radians  (6°)  is  a  reasonable  disturbance  to  the  system. 
This  does  not  mean  that  the  system  is  restricted  to  motions 
within  6°  of  the  hover  condition,  only  that  the  response 
time  for  which  we  design  applies  up  to  these  limits. 

4.  Iterative  Design  Process 

We  are  now  ready  to  form  cost  functions  for  each 
subsystem  and  determine  suitable  gain  matrices.  A  numerical 
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algorithm  for  solving  Equation  3.13  is  discussed  in  [3:83]. 
A  computer  program  for  iterating  the  solution  and  converging 
to  the  steady  state  is  described  in  [13]  and  is  used  to 
determine  the  gain  matrices  in  this  section. 

Once  a  gain  matrix  is  determined,  the  time  response 
of  the  linearized,  closed  loop  system  can  be  tested  using 
the  linear  system  computer  simulation  in  Appendix  C.  ~  The 
process  repeats  until  an  acceptable  time  response  is 
obtained. 

a.  Roll/Throttle  Subsystem 

Table  3.10  summarizes  the  results  with  various 
cost  functions.  The  values  for  Q  are  the  diagonals  of  the  Q 
matrix  which  are  the  weights  on  <p,  u,  SB  ,  5t  ,  s,  ,  and  St , 
respectively.  The  values  for  R  are  the  weights  on  ua  and  ut 
in  that  order. 

Iteration  (1)  is  with  a  Q  normalized  to  account 
for  scaling  differences  between  states.  Results  indicate 
that  the  normalized  cost  function  is  not  adequate.  A 
shorter  response  time  is  desired. 

Iterations  (2),  (3),  (4),  and  (5)  chronicle  a 

trial  and  error  design  process  where  each  succeeding  cost 
function  is  chosen  based  on  the  response  of  the  previous 
ones.  A  dominant  consideration  is  the  coupling  between  the 
throttle  and  roll  ctates.  Since  is  undesirable  to  have  the 
throttle  used  to  correct  the  roll  angle,  a  heavy  cost  is  put 
on  maintaining  the  throttle  at  the  desired  state. 
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TABLE  3.10 

Settling 

Overshoot 

Cost 

Initial 

Time 

%  (absolute) 

Deviation 

* 

Value 

U 

ID 

_ IL 

_ it 

it 

(1) 

(p  m.l 

4 

6. 

0% 

(0) 

.02 

0 

.2 

0 

(2)  _ 

.<P  =.l 

1 

.5  3. 

12% 

(0) 

.03 

-.0002 

.35 

.025  _ 

(3) 

ip  =.l 

1 

.5  2.5 

10% 

(0) 

.03 

-.004 

.35 

.06 

(4) 

<p  *.l 

1 

2. 

6% 

(0) 

.07 

-.0007  1^2 

-.007 

(5) 

<p  =.1 

1 

2. 

3% 

(-.0001)  .055  .0125 

.87 

.3 

U  =.l 

2, 

1.  (.004) 

5% 

.045  -3. 

.275 

-50 

Post  J 

functions 

> 

(1) 

(2) 

(3) 

(4) 

(5) 

Q  . 

R 

Q 

R 

Q 

R 

Q 

R 

Q 

R 

10 

[60] 

20 

M 

50 

M 

50 

fl] 

50 

r*1! 

3 

L60J 

1000 

L-iJ 

2000 

L-iJ 

2000 

L*iJ 

2000 

L*o5J 

60 

.1 

.01 

.01 

.01 

60 

.1 

.01 

.01 

.01 

30 

1 

3 

.2 

.4 

30 

.1 

.1 

.1 

.005 

Li  ^ 

-  The  final  iteration  results  in  an  acceptable 

time  response  and  satisfies  the  state  constraints.  The 
optimal  gain  matrix,  KPT  converged  in  14  iterations  at 
AT=.04  seconds  per  iteration.  Hence,  the  steady-state 
constant  gain  schedule  is  reached  after  .56  seconds.  Since 
the  settling  time  is  1  and  2  seconds,  our  assumption  that  a 
constant  gain  schedule  would  provide  the  optimal  solution  is 
valid. 

b.  Pitch/Yaw  Subsystem 

Table  3.11  lists  results  for  the  iterative 
design  of  the  pitch/yaw  subsystem  gain  matrix.  The  values 


87 


for  Q  are  the  diagonal  elements  corresponding  to  eb  ,  V>b ,  6, , 
Sr  ,  s\  ,  and  Sr  ,  respectively,  u,  and  ur  are  the  control 
inputs  in  R. 


TABLE  3.11 

Settling  Overshoot 

Cost 

Initial  Time  %  (absolute)  Deviation 

#  Value  6 _ . 

6 

*  •$,  6 -  __  $ 

Jr 

(1) 

e  ».i  2.  1.5 

10% 

(.013)  .035  -.055 

.35 

-.45 

(2) 

e  =.i  .6  l. 

8% 

(.0125)  .08  .045 

.6 

.5 

(3) 

e  =.l  .75  1.5 

3% 

(.0125)  .05  -.065 

.7 

-.7 

yf>  =.l  1.  1.5  (-.012)  2%  .065  .06 

.75 

.6 

Cost  Functions 

(1) 

(2) 

(3) 

_Q 

R 

-  % 

R 

r  % 

R 

100 

r  *1 

100 

r  *1 

150 

r  *1 

100 

LiJ 

100 

LiJ 

150 

L1] 

.01 

.01 

.01 

.01 

.01 

.01 

2 

1 

1 

2 

1 

1 

— 

The  same  procedure  as  roll/throttle  for  iterat¬ 
ing  through  various  possible  cost  functions  is  used  for  the 
pitch/yaw  subsystem.  As  in  Table  3.10,  only  selected 
iterations  are  listed.  Insight  gained  from  the  roll/ 
throttle  iterations  was  valuable  in  reducing  the  steps 
needed  for  the  pitch/yaw  design. 

The  final  gain  matrix,  KQR,  converged  in  13 
iterations,  or  .52  seconds  which  is  less  than  the  settling 


time  of  .75  seconds.  Therefore,  the  steady-state  gain 
matrix  is  an  optimal  solution  for  this  cost  function. 
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c.  Steady-State  Gains 

Table  3.12  lists  the  optimal  steady-state  gain 
matrices  resulting  from  the  iterative  design.  A  test  for 
stability  places  all  eigenvalues  of  the  discrete  time  closed 
loop  subsystems  inside  the  unit  circle.  Therefore,  the  gain 
matrix,  K,  has  succeeded  in  stabilizing  the  open  loop  plant. 


TABLE  3.12 

STEADY-STATE  GAIN  SCHEDULES 

Itoll/lhrottle 


‘  -1.99 

-.603 

-.717 

.274  1.75 

.00994 

.140 

.00016 

Kpr  “ 

-.585 

U 

111. 

-.205 

28.9  .582 

1.67 

.00962 

.137 

Pitch/Yaw 

“-1.61 

-1.61 

-.553 

.320  1.05 

-.208 

.127 

-.002641 

Kqr  * 

1.57 

-1.60 

-.305 

-.566  .107 

1.23 

.00096 

.130  i 

1 - - - i 

5.  Linear  System  Optimal  Regulator  Response 

The  graphical  results  discussed  here  are  for  the 
optimal  gains  of  Table  3.12  applied  to  the  linear  model  with 
state  constraints.  The  DSL/VS  simulation  is  listed  in 
Appendix  C . 

Figures  3.8(a),  (b) ,  (c) ,  and  (d)  are  the  regulator 
time  responses  for  initial  conditions  of  ^(0)=.l,  U(0)=.i, 
eb(0)=.l,  and  i'  '0)  =  .l,  respectively.  The  design  goals  for 
overshoot  and  settling  time  are  satisfied  for  the  linear 
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Figure  3.8  (a)  Linear  System  Regulator  Time  Response 
with  Optimal  Gains,  p(0)«.l 
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Figure  3.8  (c)  Linear  System  Regulator  Time  Response 
with  Optimal  Gains,  e ( 0) » .  1 
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Figure  3.8  (d)  Linear  System  Regulator  Time  Response 
with  Optimal  Gains,  V»(0)  =  .l 


system.  Theca  results,  for  which  the  optimal  gains  were 
designed,  are  a  reference  to  which  results  with  the  non¬ 
linear  system  can  be  compared. 


E.  THE  OPTIMAL  TRACKER 

The  optimal  gains  derived  for  the  regulator  can  be 

m.  .  «*»- 

extended  to  the  tracking  problem.  The  tracking  problem 
characterizes  the  command-response  relationship  desired  to 
steer  and  maneuver  AROO.  In  theory,  several  approaches  to 
the  tracking  problem  exist.  The  method  used  in  this  work 
utilizes  the  error  state.  The  control  law  which  applies  is 

Uc  (k)  -  -K(x(k)  -  rc  (k) )  (3.17) 

where  rc  (k)  is  the  reference  (commanded)  state  vector. 


This  control  law  is  pictured  in  Figure  3.9.  The  error  state 


— 

x(k)r - 

■={>  H  =$>  yoo 

‘ =MX)=d-L/  N>  x(k*Q=0xCk)  +  ruc(IO  \= 

Figure  3.9  Tracker  Control  Scheme 


as  it  drova  the  system  states  to  zero  in  the  regulator 
problem.  A  full  state  reference  vector  is  required  which 
means  a  commanded  value  must  be  specified  for  every  system 
state.  The  reference  vector  for  the  roll/throttle  subsystem 
becomes 

^ »^e  »^re  *^tc  »^ie  >^te  1*  ( 3  •  18][ 

and  for  pitch/yaw 

rCQR  -  [ec,tfc,QcfR e,ff.e.firc»*.ef*re]t  (3.19) 

where  the  values  of  the  reference  vectors  may  be  time 
varying  and  are  the  desired  values  of  the  final  state. 

Application  of  Eq.  3.17  to  the  linear  model  is  described 
in  Appendix  C.  Figure  3.10(a)  is  the  response  for  desired 
states  v>-e-.l.  The  response  is  within  the  desired  limita¬ 
tions  for  settling  time  (2  seconds)  and  overshoot  (10%) . 
Coupling  is  quite  pronounced  between  the  pitch  and  yaw. 

>i 

Figured.  10(b)  is  the  response  to  U«.5  and  coupling  between 
the  throttle  and  roll  heading  is  seen.  Finally,  Figure  3.10 
(c)  is  the  response  to  a  step  input  on  yaw  (tfb  =  .l).  The 
results  are  within  specification. 

F.  RESULTS  WITH  THE  NONLINEAR  SYSTEM 

The  optimal  steady-state  gains  for  the  linearized  AROD 
model  are  determined  in  the  previous  section.  The  goal  in 
designing  the  control  algorithm  is  to  provide  fast  response 
to  disturbances  and  command  inputs,  while  maintaining 
minimum  overshoot  of  the  final  state.  The  requirement  for 
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Initial  States:  vm0,  U*0,  eb-o,  ^b*0 
Desired  States:  p-.l,  U«0,  eb*.l,  Vb-0 


Figure  3.10  (a)  Linear  System  Tracker  Time  Response 
with  Optimal  Gains,  Roll  and  Pitch  Command 
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Initial  States:  <pm0,  U«0,  ©b“0,  ^-0 
Desired  States:  p«0,  U*.5,  ©b*0,  tfb*° 


Figure  3.10  (b)  Linear  System  Tracker  Time  Response 
with  Optimal  Gains,  Altitude  Rate  Command 
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Initial  States:  v o,  u-0,  eb«0,  y^-o 
Desired  States:  <p*0,  U-0,  eb«0,  ^,-.1 


I 

l 

Figure  3.10  (c)  Linear  System  Tracker  Time  Response 
with  Optimal  Gains ,  Yaw  Command 
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quickness  is  driven  by  the  erratic  behavior  of  the  full 
model  in  Chapter  II  when  it  is  disturbed  from  the  steady- 
state  pitch  or  yaw  positions.  This  section  applies  the 
linear  system  controller  to  the  nonlinear  model.  Results 
are  given  rirst  for  the  controller  applied  to  the  nonlinear 
model  without  the  pitch  and  yaw  aerodynamic  moments  and 
second  "to  the  full  model  with  all  moments  and  forces. 

1.  Nonlinear  Model  without  Aerodynamic  Moments 

Figure  3.11(a)  is  the  nonlinear  system  step  response 
for  the  heading  angle,  <p.  The  settling  time  is  1.5  seconds 
which  is  the  same  as  for  the  linear  system.  However,  the 
overshoot  is  nearly  15%  which  is  larger  than  the  linear 
case.  Figure  3.11(b)  is  the  response  to  a  desired  altitude 
rate,  U' ,  of  .1.  Note  that  for  the  nonlinear  model,  the 
altitude  rate  is  derived  from  the;  Euler  angles  and  it  is  no 
longer  approximated  by  the  body-fixed  velocity,  U.  The 
altitude  rate  settles  in  less  than  1.5  seconds,  but  the 
heading  angle,  disturbed  through  coupling,  settles  in  1.75 
seconds  with  .05  radians  overshoot.  Still,  the  system  does 
settle  quickly  and  with  zero  steady-state  error. 

Figure  3.12  is  the  system  response  for  a  desired 
state  step  input  on  yaw,  ,  of  .1.  Severe  coupling  to  <p 
and  U'  is  observed,  while  pitch,  eb ,  behaves  as  expected. 
The  reality  of  the  linearizing  assumptions  becomes  evident 
here.  The  results  simply  indicate  that  as  the  body  is 
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Initial  States:  <p*o,  U'-O,  ©b -o ,  ^-o 
Desired  States:  *-.l,  U*«0,  ^-0,  y^-o 


Figure  3.11  (a)  Nonlinear  System  Tracker  Response 
Heading  Angle  Command,  v 
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Figure  3.11  (b)  Nonlinear  System  Tracker  Response 
Altitude  Rate  Command,  U' 
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Initial  States:  *-0,  U‘«0,  ^-0,  ^-0 
Desired  States:  ¥>«0,  U*»0,  6^*0,  V»b “ •  1 


Figure  3.12  Nonlinear  System  Tracker  Response 

Yaw  Command 
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pitched  (or  yawed) ,  some  of  the  altitude  holding  thrust  will 
be  lost  and  redirected  to  a  horizontal  velocity.  Figure 
3.13  shows  that  if  the  pitch  or  yaw  is  kept  small  enough 
(here  0b«.O3  radians),  then  the  altitude  rate  can  also  be 
controlled.  Figure  3.14  goes  a  step  further  and  illustrates 
that  if  a  bounded  input  (pitching  the  body  for  1  second)  is 
applied,  then  a  speed  across  the  ground  can  be'  realized 
while  maintaining  the  desired  heading  and  altitude  rate. 

2.  Nonlinear  Model  Complete 

Results  for  the  controller  applied  to  the  full  model 
including  the  troublesome  pitch  and  yaw  aerodynamic  moments 
are  shown  here.  Figure  3.15  is  the  system  response  to  a 
step  input  on  roll  (p>.l)  A  steady-state  error  of  -.04  is 
and  is  due  to  the  downwash  swirl  moment  previously  dis¬ 
cussed.  However,  coupling  to  the  other  states  causing  the 
instability  noted  in  Chapter  II  has  been  eliminated  by  the 
controller.  Figure  3.16  is  the  result  of  a  step  input  for 
U'-*.5.  Coupling  of  the  velocity  to  the  pitch  and  yaw  axes 
causes  enough  of  a  disturbance  in  angle  of  attack  to  upset 
the  system  beyond -the  constraints  of  the  control  vanes  and 
the  controller  is  unable  to  correct  it.  To  illustrate  the 
behavior  of  the  data  to  pitching  and  yawing,  Figure  3.17 
shows  the  resulting  aerodynamic  moments,  L, ,  M„ ,  N,  , 
step  input  for  altitude  rate. 


for 


Figure  3.13  Nonlinear  System  Tracker  Response 
Yaw  Command  for  Minimum  U',  *  Deviation 
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Figure  3.14  Nonlinear  System  Tracker  Response 
Yaw  Command  is  a  2  second  step 
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Figure  3.15  Full  Nonlinear  Nodal  Tracker  Response 
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Initial  Statas:  v-Q,  U'«0,  eb «0 ,  ^  «0 
Desired  statas:  *>-0,  U'».l,  ©*-0,  ^,-0 


Figure  3.16  Full  Nonlinear  Nodal  Tracker  Response 
Altitude  Rate  Command 
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Initial  States:  ^*0,  U’-0,  eb«0,  0b«O 
Desired  States:  y>-0,  U*«.l,  ©b-0,  1>bm0 


Figure  J . 17  Aerodynamic  Moments  Response  to 
Altitude  Rate  Command 
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3.  Conclusion 


Clearly,  the  controller  is  not  adequate  to  deal  with 
such  large  discontinuities  as  the  aerodynamic  moments 
present.  The  control  vanes  are  the  limiting  factor. 
Therefore,  the  data  or  the  way  it  is  related  to  the  states 
in  the  model  is  not  correct. 

'  However,  we  concluded  in  Chapter  II  that  ~  the 
behavior  of  AROD  without  these  moments  closely  resembles  the 
gyroscopic  model  that  we  expected.  If  the  controller  is 
judged  based  on  its  performance  with  the  nonlinear  model 
without  the  pitch  and  yaw  moments,  then  the  control  design 
is  valid.  The  fact  that  the  controller  maintains  stability 
while  directing  commanded  inputs  to  steer  and  maneuver  the 
vehicle  supports  the  linear  model  approach  taken  in  this 
work. 
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IV.  RAPIQ  EBEQBBacx  M353fflft  BESISM 


A.  REQUIREMENTS  FOR  A  RF  ANTENNA 

Detailed  technical  requirements  for  the  AROD  video 
downlink  have  not  been  specified  by  the  GATERS  project 
office  as  of  this  writing.  Neither  have  the  tactical  or 
operational  requirements  been  addressed  in  detail.  However, 
GATERS  project  quarterly  review  conferences  have  undertaken 
discussion  of  this  aspect  of  AROD  and  some  consensus  has 
emerged  [2].  It  is  from  these  discussions  and  from  the 
stated  operational  mission  of  AROD  (see  Chapter  I)  that 
assumptions  are  made  on  which  to  base  the  antenna  design. 

Several  facts  are  known  about  AROD  which  pertain  to  RF 
wave  propagation. 

(1)  Cost  and  weight  are  key  factors  and  a  simple  design 
utilizing  inexpensive,  readily  available  materials 
is  most  appropriate. 

(2)  The  AROD  is  intended  to  be  man  carried  and  easily 
readied  for  launch.  Rugged,  lightweight 
construction  is  required. 

(3)  A  frequency  range  in  the  UHF  band  of  800  to  900 
megahertz  (MHz)  has  been  tentatively  identified  for 
video  signal  transmission. 

(4)  A  range  of  at  least  five  kilometers  is  required  if 
AROD  is  to  be  operated  at  the  maximum  planned 
distance  from  the  operator. 

(5)  An  omnidirectional  pattern  is  desired  for  training 
purposes  to  ensure  that  the  video  signal  can  be 
received  regardless  of  AROD's  orientation. 
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1.  * 


(6)  Angular  speed  of  the  propeller  during  hover  Is  ibout 
7200  rpm  (120  Hz)  and  a  three  blade,  carbon  fiber 
composite  propeller  is  to  be  used. 

(7)  The  body  is  made  of  a  carbon  fiber  composite  and  the 
frame  is  composed  of  aluminum.  The  forebody  is 
fiberglass. 

(8)  The  AROD  contains  an  on  board  computer,  sensors,  and 
a  radio  receiver  for  the  command  uplink. 

This  list  by  no  means  includes  all  of  the  factors  which 

<  . 

contribute  to  the  electromagnetic  character  of  AROD  nor  does 
it  address  every  possible  requirement  for  a  RF  antenna  on 
AROD.  However,  these  requirements  allow  a  design  goal  to  be 
determined  based  on  the  following  assumptions: 

(1)  The  120  Hz  disturbance  from  the  rotating  propeller 
and  the  command  uplink  RF  signal  will  not  interfere 
with  the  800  to  900  MHz  transmission  signal.  An 
order  of  magnitude  between  signals  is  a  common  rule 
for  non-interfering  signals.  The  propeller  (7200 
rpmj  and  the  command  uplink  (72  MHz)  meet  this 
constraint. 

(2)  For  the  purpose  of  electromagnetic  analysis,  the 
AROD  can  be  modelled  as  the  cylinder  pictured  in 

_ Figure  4.1.  The  cylinder  has  the  outside 

dimensions  of  the  AROD  main  body  surface  shown  in 
Figure  1.1. 

(3)  The  carbon  fiber  body  is  highly  conductive.  AROD 
project  engineers  have  stated  that  the  composite 
used  has  high  conductivity  characteristics. 

(4)  The  carbon  fiber  body  will  essentially  shield  out 
electromagnetic  interference  from  the  on  board 
computer  and  sensors. 

(5)  The  antenna  will  conform  to  the  body  of  AROD  as  much 
as  possible  to  prevent  breakage . while  being  carried 
and  snagging  during  operation  in  closed-in  areas. 

(6)  The  antenna  will  be  fixed  to  facilitate  simple 
launch  preparation. 
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B.  DESIGN  SELECTION 

A  -logical  point  to  start  tha  choice  of  antenna  is  with 
devices  already  in  use  on  similar  vehicles.  AROD  is  most 
similar  to  a  satellite  in  that  it  is  free  flying  yet  not 
fast  moving  and  aerodynamic  influences  are  negligible.  One 
UHF  satellite  antenna  in  common  use  is  the  mesh  array 
commonly  referred  to  as  the  umbrella  [14].  However,  the 
umbrella  antenna  is  neither  rugged  nor  omnidirectional  and 
would  be  a  poor  choice  for  AROD. 

Stripline  antennas  are  arrays  of  slotted  strips  which 
are  surface  mounted  to  a  vehicle.  The  slotted  strips  are 
cut  for  specific  frequencies  and  are  generally  not  tuneable. 
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These  antennas  ara  lightweight,  abla  to  ba  flush  mounted  and 
could  ba  aada  ruggad  and  durable  for  AROD.  Research  for 
this  report  has  found  only  L-band  and  S-band  variations  of 
the  stripline  satellite  antenna  in  use.  However,  the 
concept  of  a  flush  mounted  stripline  can  potentially  be 
exploited  in  the  development  of  the  UHF  antenna  for  AROD. 
[15]  A  simple  antenna  to  employ  is  a  dipole  extending 
from  the  top  of  AROD  or  out  the  bottom.  At  a  transmission 
frequency  of  800  MHz,  the  wavelength  is 


14.76  inches 


(4.1) 


where  c  -  1.18  x  101 0  inches/second. 


A  half-wave  dipole,  for  example  would  extend  7.38  inches 
above  or  below  AROD.  This  design  is  not  considered  here. 
Instead,  an  antenna  which  conforms  to  the  AROD  body  is 
developed  to  satisfy  two  of  the  needs  listed  above:  the 
need  for  a  rugged,  fixed  antenna  and  the  presence  of  a 
conducting  surface  (the  body) . 

The  need  for  a  rugged,  fixed  antenna  is  obvious  and  will 
not  be  discussed  further.  The  conductivity  of  the  AROD 
body,  however,  can  possibly  enhance  the  conformal  antenna's 
radiation  pattern  characteristics  if  the  body  behaves  like  a 
reflecting  ground  plane.  This  potential  and  basic  antenna's 
which  can  take  advantage  of  the  reflecting  body  are  dis¬ 
cussed  in  the  next  section. 


C.  THEORETICAL  DEVELOPMENT 

Two  basic  designs  for  the  UHF  antenna  are  considered 
here.  Referring  to  Figure  4.1,  the  first  is  a  vertically 
polarized  dipole  or  dipole  array  and  ,  the  second  is  a 
horizontally  polarized  loop  antenna.  The  dipole  array  and 
the  loop  antenna  are  both  very  well  developed  in  References 
16  and'  17 .  This  work  does  not  endeavor  to  recreate  or 
expand  on  oither  books'  theoretical  development  and  the 
reader  is  encouraged  to  consult  these  sources  if  such 
information  is  desired.  Instead,  a  brief  explanation  of  the 
important  characteristics  of  the  dipole  and  the  loop  antenna 
is  offered  as  a  preface  to  the  practical  considerations 
which  drive  the  design  process. 

1.  Dipole  Array 

a.  The  Half-Wave  Dipole 

An  omnidirectional  far  field  pattern  is  the 
goal.  For  this  reason,  a  dipole  which  radiates  an  omni¬ 
directional  pattern  is  desired.  An  obvious  choice  is  the 
half-wave  dipole.  Fed  at  the  center,  the  current  on  a  half¬ 
wave  dipole  will  decrease  sinusoidally  as  it  nears  the  end 
points. 

A  single  half-wave  dipole  yields  the  far  field 
pattern  of  Figure  4.2.  The  horizontal  (azimuth)  pattern  is 
omnidirectional.  If  the  X  direction  of  the  dipole  cor¬ 
responds  to  the  X  direction  of  AROD  then  it  is  the  horizon¬ 
tal  pattern  of  the  dipole  which  is  of  greatest  interest  in 
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Figure  4.2  Far  Fiald  Patterns  of  a  Half-Wave  Dipole 

obtaining  an  acceptable  solution  for  AROD.  Therefore,  the 
emphasis  is  placed  on  obtaining  an  acceptable  horizontal  far 
field  pattern  from  subsequent  antenna  designs. 

b.  Circular  Array  of  Half-Wave  Dipoles 

A  dipole  which  is  close  to  the  AROD  body  will 
have  radiation  characteristics  far  different  from  the  free 
space  or  free  standing  half-wave  dipole.  The  conducting 
body  will,  at  the  very  least,  block  the  dipole  pattern  in  a 
direction  "behind"  or  toward  the  body  from  the  dipole.  An 
array  of  dipoles  which  conforms  to  and  encircles  the  body 


will  allow  a  signal  to  radiate  through  360  dagraas  in 
azimuth  (*>) .  Such  an  array  is  callad  a  circular  or  ring 
array  [16:p.  352].  A  circular  array  of  identical,  aqually 
spacad  half-wava  dipolao  aach  axcitad  at  its  cantar  will 
rasult  in  a  non-diractional  far  fiald  pattam.  Figura  4.3 
is  tha  rasulting  far  fiald  gain  pattam  of  a  four  dipola 
ring  array.  Although  tha  pattam  is  plaguad  by  peaks  and 
valleys  craating  a  wida  dynamic  ranga  in  diractiva  gain,  it 
is  a  good  starting  point  for  tha  dasign. 

The  theory  and  analysis  for  circular  arrays  is 
well  developed  and  an  array  yielding  an  acceptable  omni¬ 
directional  pattern  can  be  derived  using  the  procedure  given 
in  [I6:p.  354].  However,  no  procedure  is  found  for  an  array 
with  an  obstructing  body  in  its  center.  For  this  reason,  a 
model  of  the  AROD  body  which  approximates  its  electromag¬ 
netic  characteristics  is  useful.  To  such  a  model,  'various 
arrays  can  be  added  and  tha  rasulting  far  field  patterns  can 
be  analyzed  for  adequacy. 

2.  Loop  Antannaa 

The  loop  antenna  is  an  attractive  alternative  for 
the  AROD  application.  A  loop  YKUld  naturally  conform  to  the 
body  circumference  and  could  be  constructed  easily.  Loop 
antennas  are  thoroughly  discussed  in  Reference  17.  However, 
analytical  results  are  readily  derived  only  for  ’'small"  or 
"moderate"  sized  loops.  A  moderate  size  loop  is  defined  as 

^  b  <  2.5  (4.2) 
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Figure  *i.3  Horizontal  Far  Field  Pattern  -  4  Dipole  Ring 
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where 

Pi  -  3.1415962, 
b  is  ths  radius  of  ths  loop, 

A  is  ths  wavslsngth. 

A  circumferential  loop  on  Ai'OD  has  a  radius  b  £  15  inchas. 
At  a  transmission  frsqusncy  o.:  800  KHz,  ths  wavslsngth  A  « 

14.76  inches  and  aquation  4.2  bscomss: 

*  .  ,  ^ 

b  >  6.39  (f  -  800  MHz)  (4.3) 

Clsarly,  AROD  doss  not  satisfy  ths  requirements  for  a 
moderate  size  loop. 

The  difficulty  with  deriving  analytical  results  for 
larger  than  moderate  size  loops  lies  in  evaluating  the 
Maxwell's  Equations.  Larger  loops  such  as  an  antenna  around 
the  AROD  body  require  a  much  mors  extensive  mathematical 
development  to  evaluate  ths  equations  than  is  desired  here. 
Therefore,  a  numerical  technique  with  the  aid  of  a  computer 
is  best  suited  for  evaluating  the  radiation  pattern  of  the 
proposed  design. 

Figure  4.4  is  the  far  field  gain  pattern  of  a  single 
loop  antenna  with  a  radius  b  ■  15  inches.  While  this 
pattern  is  apparently  directional  and  inadequate  for  the 
AROD  application,  combined  with  the  AROD  body  model,  a 
better  pattern  may  result. 

3.  sraunfl  Elans  BE tests 

Perfect  ground  planes  are  defined  as  infinite, 
plane,  and  perfect  conductors.  Practical  approximations  to 
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Figure  4.4  Horizontal  Far  Field  Pattern-  Single  Loop 


a  perfect  ground  plana  ara  planas  much  largar  than  1*  and 
good  conductors  on  tha  ordar  of  10?  mhos/meter.  [16 ip.  239] 
Claarly,  a  cylindrical  surfaca  only  slightly  largar  than  2A 
at  its  graatast  dimansion  can  hardly  be  considarad  infinita 
and  plana.  Navarthalass,  tha  fact  that  tha  AROD  body  doas 
conduct  suggests  that  soma  link  to  ground  plana  theory  can 
ba  achieved. 

Ground  plana  theory  is  tha  study  of  reflections.  An 
ideal  dipole  oriented  parallel  to  a  perfect  ground  plana  has 
an  image  in  tha  plana  which  is  directed  out  from  tha  plana. 
Tha  affect  is  to  double  tha  image  radiated  outward  and 
perpendicular  to  tha  plana.  In  tha  case  of  AROD,  tha  image 
is  increased  outward  from  tha  body  in  tha  horisontal  plana. 
Further,  since  tha  body  surfaca  is  curved,  tha  image  will  ba 
reflected  at  increasingly  obtuse  angles  as  distance  from  tha 
dipole  increases.'  Tha  object  of  modeling  tha  AROD  body  is 
to  determine  if  this  pattern  of  reflection  can  ba  used  to 
enhance  the  far  field  patterns  of  Figure  4.3  and  4.4  and 
obtain  an  adequate  antenna  design. 

a.  Ground  Planes  and  Phase 

Before  proceeding,  an  important  consideration  is 
phasing  of  the  signal.  The  800  to  900  MHz  waveform  planned 
for  the  video  signal  carries  with  it  a  period  and  phase. 
The  phase  of  this  waveform  will  vary  as  the  distance  from 
the  antenna  increases.  Reflection  theory  stipulates  that 
the  distance  a  waveform  travels  to  and  from  a  reflecting 
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surface  can  ba  viewed  as  a  straight  lina  distanca  front  the 
original  source,  in  this  casa  the  antanna.  Thus,  a  phase 
change  is  associatad  with  any  reflection.  This  phasa  changa 
can  act  to  complicate  or,  at  180  dagraas,  it  can  cancal  tha 
original  signal.  To  minimise  any  such  complications,  tha 
antanna  should  ba  positionad  flush  to  tha  raflacting  surfaca 
or  at  'a  distanca  rasulting  in  zero  phasa  changa  in  tha 
raflactad  signal.  Furthar  work,  not  considarad  hara,  night 
utilize  phasing  as  a  naans  to  staar  tha  pattara  as  dasirad. 
b.  Ground  Pl&nas  and  tha  Loop  Antanna 

Ona  vary  important  rasult  of  uodalling  tha  body 
as  a  conductor  and  approxinating  it  as  a  ground  plana  is 
that  it  will  significantly  altar  tha  charactaristics  of  tha 
loop  antanna.  In  fact,  a  circumfarantial  loop  around  tha 
cylindrical  body  night  altaroativaly  ba  considarad  a  lina 
sourca  ovar  a  ground  plana.  Tha  "lina"  is  the  loop  which 
carrias  a  currant  from  its  sourca.  Unlike  a  lina  sourca, 
however,  this  currant  would  ba  seen  on  tha  loop  from  two 
different  directions  since  tha  "line"  ends  are  connected 
together,  -with  these  important  differences  in  mind,  line 
source  effects  will  play  a  major  part  in  shaping  the  numer¬ 
ically  generated  far  field  patterns  in  the  next  section. 

D.  BODY  MODELLING  AND  COMPUTER  AIDED  DESIGN  RESULTS 

The  complexity  of  the  body  conforming  antenna  makes  use 
of  a  computer  algorithm  to  model  it  very  practical. 
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Although  a  simple  array  of  dipoles  lends  Itself  to  analyti¬ 
cal  solution,  a  large  loop  the  size  of  the  one  proposed  does 
not.  A  computer  program  written  at  the  Naval  Ocean  Systems 
Center  entitled  Numerical  Electromagnetics  Code  (NEC)  is 
available  on  the  Naval  Postgraduate  School's  mainframe 
computer  [18].  This  code  is  used  to  obtain  the  far  field 
patterns  pictured  in  this  work  and  is  very  suitable  for  this 
study . 

The  AROD  body,  although  it  can  be  pictured  as  a  simple, 
conducting  cylinder  (Figure  4.1),  is  a  very  complex  electro¬ 
magnetic  structure.  One  approach  to  simulate  such,  a 
cylinder  is  to  construct  a  wire  grid  outline.  Another 
method  is  to  utilize  a  continuous  surface  model.  Both 
techniques  are  discussed  in  [16]  and  both  can  he  modelled 
with  NEC.  In  fact,  because  of  the  symmetry  of  the  AROD  body 
and  the  use  of  symmetry  possible  in  NEC,  the  continuous 
surface  model  is  somewhat  less  cumbersome  to  create.  The 
modelling  process  and  the  resulting  computer  code  are  the 
subject  of  Appendix  D. 

-  1-  Results  with  Dipole  Arrays 

The  dipole  arrays  previously  discussed  are  simulated 
here  with  the  AROD  body  model.  While  array  effects  are 
expected  to  influence  the  resulting  far  field  pattern,  the 
conducting  surface  of  the  AROD  will  contribute  significantly 
as  well.  All  runs  utilize  a  single  center  feed  on  a  .1  inch 


diameter  wire  which  is  divided  into  11  equal  segments  to 


make  a  half-wave  dipole.  The  first  antenna  simulated  with 
the  body  model  is  the  four  dipole  ring  array  pictured  in 
Figure  4.5.  Note  that  the  array  radius  b  »  15.5  inches. 


Figure  4.5  Body  Model  with  4  Dipole  Ring  Array 

The  half  inch  beyond  the  body  radius  of  15  inches  is  to 
account  for  practical  considerations  of  mounting  the  dipoles 
and  to  prevent  "electrical  contact"  with  the  body  in  the 
computer  model.  Results  of  the  four  dipole  simulation  are 
shown  in  Figure  4 . 6  and  should  be  compared  with  the  free 
standing  ring  array  of  Figure  4.3.  While  the  resulting 
pattern  is  hardly  adequate,  it  does  demonstrate  the  benefi¬ 
cial  effects  of  the  reflecting  AROD  body.  Directivity  has 
been  reduced . 


Figure  4.6  indicates  that  the  body  has  effectively 
duplicated  the  pattern  of  each  dipole  and  doubled  the 
radiating  horizontal  pattern.  Since  the  pattern  multiplica¬ 
tion  depends  largely  on  the  distance  between  radiating 
sources,  an  analysis  of  the  spacing  is  revealing.  The 
circumference  of  the  30  inch  diameter  AROD  is  94.2  inches. 
The  antenna  radius  b  *  15.5  inches  yields  an  array  radius  of 
97.38  inches.  The  four  dipoles,  then,  are  each  24.3  inches 
apart  along  the  circumference.  In  wavelengths  this  distance 
is  1.65A.  This  is  not  the  optimal  spacing.  The  spacing 
sought  is  pictured  in  Figure  4.7.  If  the  cylindrical 


Figure  4.7  Pattern  Multiplication  of  Dipoles 
Above  a  Ground  Plane 

surface  can  be  ignored  and  instead  thought  of  as  a  plane, 
then  the  interaction  of  the  dipole  patterns  can  be  supposed. 
Figure  4.2  showed  the  horizontal  pattern  resulting  from  a 
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single  half  wave  dipole.  If  the  surface  reflects,  then  the 
patterns  from  two  adjacent  dipoles  will  combine.  Dipoles  a 
half  wavelength  apart  will  result  in  patterns  adding  exactly 
in  phase  producing  sharp  peaks  and  valleys  as  seen  in  Figure 
4.7.  Patterns  from  dipoles  a  quarter  wavelength  apart  will 
tend  to  smooth  each  other.  Therefore,  spacing  close  to  .251 
is  sought  to  reduce  the  dynamic  range.  Equal  spacing  is 
needed  to  maintain  a  non-directional  array  as  previously 
discussed. 

A  circumference  of  97.38  inches  and  wavelength  A  = 
14.76  inches  yields  no  whole  multiple  of  .25.  TablU?  4.1 
lists  the  multiples  of  quarter  wavelength  spacing 


TABLE  4.1 

Approximate 

Number  of 

Actual 

— 

Spacing 

Dipoles 

Spacing 

A/4 

26 

3.7  inches 

3  A/ 4 

9 

10.8  inches 

5A/4 

5 

19 . 3  inches 

7A/4 

4 

24.3  inches 

circumference-  97.8 

inches 

A/4- 

3.69  inches  (800  Mhz) 

possible  on  the  body.  Four  dipoles  did  not  yield  an 
acceptable  pattern.  Figure  4.8  is  the  result  of  five 
equally  spaced  half-wave  dipoles.  The  pattern  is  clearly 


126 


more  smoothed  although  more  complex  than  the  four  dipole 
case.  The  horizontal,  vertical  and  total  patterns  are  shown 
in  Figure  4.8,  but  the  total  pattern  (solid  line)  is  of  most 
interest  here.  In  this  case  the  smoothing  effects  of  the 
reflecting  body  has  reduced  significantly  the  sharp  drop  off 
in  directive  gain  characterized  by  the  four  dipole  example. 

*  The  case  of  a  nine  dipole  array  deserves  some 
practical  discussion  about  feed  points.  The  problem  of 
feeding  multiple  points  in  phase  from  one  transmitter,  as  in 
the  case  of  AROD,  is  easy  to  overcome  on  the  computer  but 
more  difficult  with  actual  components.  Finite  distances 
exist  between  transmitter  and  feed  point  and  phasing  is 
simplified  with  fewer  feed  points.  Additionally,  power 
losses  associated  with  resistive  antennas  may  become  too 
great  with  multiple  loads.  Therefore,  for  the  nine  dipole 
array,  the  simulation  will  consider  only  three  equally 
spaced  feed  points. 

Figure  4.9  shows  the  results  of  this  arrangement. 
While  not  a  very  omnidirectional  pattern,  the  nine  dipole 
case  reveals  that  reducing  the  number  of  feed  points  does 
not  increase  the  directivity  and, to  the  contrary,  reduces 
the  severity  of  the  peak  to  valley  gain  difference  of 
dynamic  range.  With  this  observation  in  mind,  the  five 
dipole  antenna  which  produced  a  smoother  pattern  may  be 
enhanced  with  fewer  feed  points.  Figure  4.10(a)  is  the 
result  with  only  one  feed  point  at  ■  0  on  the  five  dipole 
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array.  Tha  far  fiald  pattern  ia  greatly  improved  over  the 
five  feed  case  and  the  dynamic  range  is  greatly  reduced. 
Figures  4.10(b)  and  (c)  are  vertical  cuts  over  zenith  (6) 
for  two  different  azimuth  angles.  The  patterns  are  non- 
directional  and  show  some  promise  as  a  suitable  antenna. 

2.  Results  with  Loop  Antennas 

'  Two  variants  of  the  loop  antenna  with  AROD  body 
model  are  studied  here.  As  in  the  case  of  the  dipole 
arrays,  reflections  play  a  major  role  in  improving  the  far 
field  pattern  over  the  free  space  results  of  Figure  4.4. 

The  loop  is  a  single,  circumferential  wire  pictured 
in  Figure  4.11.  The  first  variant  has  one  feed  point  at  <p  - 
0.  The  horizontal  pattern  which  results  is  shown  in  Figure 
4.12(a)  and  differs  markedly  from  the  free  standing  loop 
pattern.  While  the  gain  peaks  of  the  free  standing  pattern 
can  be  detected  in  Figure  4.12(a)  at  <p  -  0,  60,  120,  240  and 
300  degrees,  the  other  peaks  which  give  a  desirable  non- 
directional  far  field  are  not  a  result  of  the  loop  theory 
analysis.  However,  if  the  antenna  is  viewed  as  a  line 
source  as  previously  discussed,  this  pattern  is  essentially 
an  800  MHz  wave  propagating  along  the  wire.  A  peak  is  seen 
every  quarter  wavelength  with  only  slight  variation  at  <p  » 
180°  probably  due  to  opposing  currents  meeting  at  this 
point.  The  currents  are  not  in  phase  when  they  meet  since 
the  distance  travelled  (half  the  circumference)  is  not  a 
whole  multiple  of  .  5A.  The  97.38  inch  circumference  fits 
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Figure  4.10  (a)  Horizontal  Far  Field 
5  Dipoles  with  Body 


Figure  4.10  (b)  Vertical  Far  Field  (p-0) 
5  Dipoles  with  Body 
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Figure  4.10  (c)  Vertical  Far  Field  (*-45° ) 
5  Dipoles  with  Body 
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Figure  4.11  Single  Loop  and  Body  Nodal 
a  6A  line  source  with  a  feed  at  p  •  0.  The  peaks  which  are 
seen  on  the  far  field  pattern  of  Figure  4.12(a)  can  be 
viewed- as  the  peaks  which  result  from  the  line  source. 

The  far  field  pattern  of  the  single  feed  loop  is 
further  described  in  Figures  4.12(b)  and  (c) .  For  the  case 
when  ip  »  0,  the  zenith  or  vertical  pattern  resembles  that  of 
a  long  (6A)  lines  source  endfire  array  as  described  in 
[16:p.  180].  Discrepancies  between  the  theoretical  endfire 
far  field  and  that  pictured  here  are  the  result  of  (1)  the 
fact  that  this  antenna  is  a  loop,  (2)  an  actual  line  source 
lenght  of  6.5A,  and  (3)  opposing  currents  due  to  the  shape 
of  the  wire  as  a  loop. 
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Figure  4.12  (a)  Horizontal  Far  Field  -  1  Loop  1  Feed 


4.12  (b)  Vertical  Far  Field  (<p-0) 
One  Loop  -  One  Feed 


Figure  4.12  (c)  Vertical  Far  Field  (^=45) 
One  Loop  -  One  Feed 
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The  second  variant  is  a  single  loop  with  two  feed 
points  fed  out  of  phase.  The  feed  points  are  at  and 
¥>=180  degrees.  The  horizontal  pattern  is  shown  in  Figure 
4.13  and  has  more  directivity  than  the  single  feed  case. 

E.  ANALYSIS  &  DESIGN  CONCLUSIONS 
1.  ‘  Analysis 

The  goal  of  the  design  process  was  to  model  the  AROD 
body  and  determine  some  basic  antennas  that  yield  omnidirec¬ 
tional  horizontal  far  field  patterns.  An  iterative  approach 
utilizing  numerical  techniques  in  the  form  of  NEC  was  used 
to  screen  various  designs.  Only  a  few  of  these  designs  were 
discussed. 

The  physical  requirement  most  restrictive  to  the 
design  was  that  it  conform  to  the  cylindrical  AROD  body. 
This  constraint  led  to  the  choice  of  half-wave  dipole  arrays 
and  a  loop.  The  study  did  not  result  in  a  perfectly 
omnidirectional  antenna ,  but  two  designs  which  show  poten¬ 
tial  for  refinement  may  have  been  found.  Adequacy  in  this 
case  is  judged  by  the  dynamic  range  that  a  radio  receiver  is 
required  to  have  in  order  to  sense  an  800  MHz  signal 
regardless  of  che  orientation,  in  azimuth,  of  AROD.  In 
other  words,  the  horizontal  far  field  pattern  of  the  AROD 
must  have  a  minimum  dynamic  excursion  between  peaks. 


Figure  4.13  Horizontal  Far  Field  (0=0) 
One  Loop  -  Two  Feed 
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The  two  designs  which  present  the  best  potential  in 
for  refinement  are  the  single  loop,  single  feed  and  the  five 
dipole,  single  feed  array.  This  conclusion  is  a  result  of 
their  overall,  relatively  smooth  horizontal  pattern.  Both 
of  these  antennas  are  described  in  terms  of  horizontal  and 
vertical  far  field  patterns  in  Figures  4.12  and  4.10. 

The  most  important  assumption  which  was  made, 
relative  to  the  resulting  patterns,  is  that  the  AROD  body  is 
a  good  conductor.  This  assumption  is  based  on  the  electri¬ 
cal  characteristics  of  the  carbon  fiber  composite  used  to 
construct  the  body.  It  should  be  noted  that  this  conducting 
shield  is  needed  to  provide  the  reflecting  surface  exploited 
by  the  antenna  designs  as  well  as  to  shield  frequency 
interference  between  the  antenna  and  the  electronics 
internal  to  AROD. 


v.  SSHCIflSIgHS 


A.  DYNAMIC  MODEL 

A  computer  model  of  the  ducted  fan  hovering  vehicle, 
AROD,  was  developed  by  considering  the  vehicle  as  both  a 
gyroscope  and  as  a  flying  vehicle.  Experimental  data  was 
integrated  with  the  model  to  simulate  the  forces  and  moments 
which  act  on  AROD  in  flight.  Problems  with  the  application 
of  some  of  the  data  to  the  computer  model  were  experienced. 
However,  the  model  without  the  problematic  data  behaved  as 
expected  of  such  a  device. 


B.  CONTROL 

A  linearized  approximation  of  the  nonlinear  model  was 
developed  so  that  optimal  control  could  be  used  to  obtain  a 
steady-state  gain  matrix  with  the  ultimate  goal  of  controll¬ 
ing  the  nonlinear  computer  model.  Optimal  control  proved  a 
viable  means  of  determining  a  constant  gain  schedule  for  the 
multiple- input,  multiple-output  linear  system.  The  result¬ 
ing  controller  also  found  application  to  the  nonlinear  model 
with  the  problematic  data  removed. 

C.  ANTENNA 

Two  basic  antenna  designs  were  analyzed  for  use  as  an 
omnidirectional  antenna  for  a  UHF  video  downlink  on  the 
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AROD.  A  five  dipole  array  and  a  single  loop  antenna  were 
found  to  have  the  best  potential  of  the  designs  which  were 
tested  using  a  computer  simulation.  The  body  of  the  AROD 
was  found  to  significantly  alter  the  directivity  charac¬ 
teristics  of  the  antennas. 


D.  FURTHER  WORK 

In  the  areas  of  dynamic  modelling  and  control  of  AROD: 

(1)  The  aerodynamic  data  and  its  application  should  be 
verified. 

(2)  A  training  simulator  based  on  the  linear  model  can 
be  implemented  utilizing  the  gain  schedule  derived 
here  on  small,  readily  available  computers. 

In  the  area  of  a  UHF  antenna: 

(1)  Modelling  the  entire  AROD  including  the  forebody  and 
testing  the  designs  for  variation  from  these 
results . 

(2)  Modelling  multiple  loop  antennas  with  differing 
phasing  in  an  effort  to  steer  the  vertical  pattern. 

(3)  Analyze  the  antennas  for  power  requirements  as 
weight  of  batteries  and  power  supplies  is  critical 
to  the  AROD. 

(4)  Test  the  designs  on  the  AROD  prototype  for 
suitability. 
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APPENDIX  A 


EXPERIMENTAL  DATA 


Empirical  data  received  from  Sandia  National 
Laboratories  is  listed  here.  This  data  was  obtained 
experimentally  through  measurements  and  wind  tunnel  testing. 


1.  Weights  and  Moments  of  Inertia 


Description 

Name 

Weight  of  Prototype 

WtPro 

Weight  of  Model 

WtMod 

Moments  of  Inertia: 

about  X-axis 

Ix 

about  Y-axis 

ly 

about  Z-axis 

value (units) 
85  (lbs) 

76.5  (lbs) 


1.5246  (ft2  lbs) 
1.6767  (ft2  lbs) 
1.6684  (ft"  lbs) 


Products  of  Inertia:  (assumed)  1^- 

Moment  of  Inertia  of  Propeller 
about  X-axis  Ir 

about  Y-axis (assumed) 
about  Z-axis (assumed) 


.015012  (ft2  lbs) 

0 

0 


Propeller  Weight (derived)  PropWt  4.50017  (lbs) 


2.  Hardware  _ Characteristics 

Propeller  Efficiency 
Minimum  Thrust  of  Engine 
Maximum  Thrust  of  Engine 
Maximum  Throttle  Setting 
Maximum  Change  in  Throttle 
Propeller  Speed  at  Hover 
Propeller  Compensation 
Max  Deflection  of  Vanes 
Max  Def  Rate  of  Vanes 


Pref 

10.472 

(rad/ lb-sec) 

Thrmin 

35 

(lbs) 

Thrmax 

115 

(lbs) 

TAmax 

RTmax 

100 

100 

(lbs/sec) 
(lbs/ sec5  ) 

HovRat 

7200 

(rpm) 

RHo 

.0019 

(ft- lb*) 

MaxDef 

.5236 

(rad) 

MaxRatDef 

.87266 

(rad/s^c) 
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Daacr lotion 


Vttlutfunif ) 


Control  Surface  Effectiveness 
Aileron 
Ruddar 
Elevator 


N 


;aaff 

eeff 

raff 


-21.29  (aac  J) 
-14.51  (aac“f ) 
-16.68  (aac  *) 


3.  Aai 


iaais. 


Tha  data  liatad  in  tha  tabla  balow  ia  obtainad  from  wind 


tunnal  tests .  Tha  manner  in  which  it  is  appliad  in  the' 
simulation  (relationship  to  forcaa  and  moments  on  tha  AROD 
body)  is  given  in  tha  next  section. 

Angle  of  Attack,  AOATOT  -  +  fi*  ) 


ttime  55 _ 

_ 55 

55 

55 

25 _ 

_ 25  55  55 

CRdel  — 

— 

.0807 

.0743 

.0712  .881  .085 

Rslope  .1514 

.1533 

.1673 

.1977 

.1834 

.1764  .2041  .20 

Req  9.293 

8.82 

8.115 

7.848 

7.347 

7.134  5.011  0.0 

CPdel  — 

— 

— 

.3106 

.3081 

.3305  .3515  0.0 

Pslope  .5533 

.512.. 

.6974 

.7612 

.7603 

.8189  .8853  .90 

Peq  43.14 

31.00 

29.58 

26.74 

24.55 

20.43  15.17  0.0 

CYdel  — 

— 

.0024 

-.002 

.0090  .0120  0.0 

Yslope  .0609 

.0663 

.  0215 

.0060 

-.004 

.0223  .0302  .04 

Yeq  . 555 

.576 

.006 

-.09 

-1.63 

-1.77  -1.74  -1.8 

VanEff  — 

— — 

1.5 

1.5 

1.5  1.9  1.0 

Veq  51.62 

44.96 

38.56 

32.94 

27.76 

21.83  13.31  0.0 

CLdel  — 

— 

.2428 

.2391 

.1907  .1734  .1046 

Lslope  . 7954 

.7183 

.6598 

.5949 

.5901 

.4726  .4366  0.0 

Leq  85. 

85. 

85. 

85. 

85. 

85.  85.  85. 

CDdel  — 

— 

—— 

.5754 

.6052 

.5922  .4966.  3890 

Dslope  1.302 

1.313 

1.371 

1.410 

1.493 

1.467  1.251  1.2 

Oeq  0.0 

0.0 

0.0 

0.0 

0.0 

0.0  0.0  0.0 

CSdel 

.0091 

.0258 

.0331  .0318  0.0 

CSslopa  .0801  .0801 

.0801 

.0801 

.0801 

.0801  .0801  .0801 

Seq  0.0 

0.0 

0.0 

0.0 

0.0 

0.0  0.0  0.0 

NOTE:  ( — )  indicates  that  no  data  under  this  heading  was 
received  for  this  angle. 
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4.  Aerodynamic  Fore—  and  Moments 

The  aerodynamic  data  listed  in  paragraph  3 ,  above, 
results  in  forces  moments  about  the  velocity  axes,  Vtot. 

The  relationships  developed  by  the  engineers  who  obtained 
the  data  are  reproduced  here.  Some  convenient  groupings  of 
constants  and  parameters  of  the  system  are  given  here. 

'  WtRat  - 

,rcRa*  wtPro 


Vdelta  -  Vtot  -  v'  (WtRat)  *Vaq (AOATOT) 

RHovA  -  RHo*Pi*  (WtRat) *Veq( AOATOT) 
wr  ■  angular  speed  of  propeller  (rad/sec) 
UseOel  »  -«r*RHovA*WtRat 
UseEq  -  (1  +  Jr*wr)*WtRat 

UseS Ip  -  WtRat* (VtQt  -  (WtRat) *Vaq (AOATOT) 
The  Forces  and  Momenta  which  result  from  the 
aerodynamic  data  are  given  here. 

FORCES  Jwith  respect  to  Vtot* 

Puerto  Body  Aerodynamics 

Fai«  UseDel*CLdel (AOATOT)  +  UseSlp*Lslope (AOATOT) 
Fad*  UseDel*CDdel (AOATOT)  +  UseSlp*Dslope (AOATOT) 
Fas"  UseOel *CSdel(AOATOA)  +  UseSlp*Sslope 
Due_to  Thrust  Force  Aerodynamics 
Ftl*  UseEq*Leq 

Total  Aerodynamic  Forces 
Fl"  Fal  +  Ftl 
Fd*  Fad 
Fs  *  Fas 
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Rar-  Ua*D*l*CRd*l (AOATOT)  +  Ua#Slp*Ralop* (AOATOT) 
Pap-  Ua*Dal*CPd*l (AOATOT)  +  Ua*Slp* Pal op* (AOATOT) 
Yay"  Ua*D«l*CYd*l (AOATOA)  +  Ua*Slp*Yalop* (AOATOT) 


•?r 


^vrw» 


ittr-  Ua*Eq*R*q (AOATOT) 
Ptp-  Ua*Eq*Paq (AOATOT) 
Yty-  Ua*Eq*Y*q (AOATOT) 


Rr“  Rar  +  Rtr 
Pp“  P®p  +  Ptp 
Yy«  Yay  +  Yty 


148 


The  following  source  ooda  is  a  Dynamic  Simulation  Language  (DSL) 
program  to  simulate  the  APOD  aa  a  gyroaoopa.  This  simulation  is 
discussed  in  Chaptar  H. 


TITLE  BASIC  AROD  GYRO  SIMULATION 

CONST  IRX-69.52, I  XX- 7063. 39, l Y Y-7768 . 2 2 , l 2Z -7729 . 58 , URX -784 . 
• 

INIT 

P0-  0 . 

00-0 . 

80-  0. 

MX  •  0. 

MZ  -  0. 

* 

DYNAMIC 

IF  (TINE  .LT.  TO)  THEN 
MY  ■  0 
(L4f 

***  STEP  INPUT  *** 

T1  -  TO  ♦  OUR 

MY  -  AMP*(STEP(TO)  -  STEP(TI)) 

END  IF 

* 

DERIV 

HX  «  I RX*URX  ♦  IXX*P 
NY  «  I Y Y*Q 
HZ  •  I ZZ*R 

PO  •(•1.*HZ/IXX)*Q  ♦  (HY/IXX)*R  «•  MX/IXX 
00  «(HZ/IYY)*P  •  <HX/IYY)*R  ♦  MY/ 1 Y Y 
RO  «<  •  1  .*HY/IZZ)*P  ♦  ( HX/I ZZ)*0  ♦  M2/IZZ 
P  -  INTGRL(PO,PD) 

0  >  INTGRL(QO.QD) 

R  -  INTGRL(RO.RO) 

PARAN  TO- 1  ,  OUR-1 ,  AMP-10000 
CONTROL  FINTIM-5,  OELT-.OS 
SAVE  .05,  P , Q , R , M Y 
PRINT  .1,MY,P,0,R 
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ili  ■  v  *  ■_  w*  u  * 


GRAPH  (G0,DE-TEK618)  T  I  ME ( N I -5 , UN«S EC )  MY ( UN- 1  LB - l N ( SQ  )  •  ) 

LABEL  (GO)  INPUT  TORQUE 

GRAPH  (G1 >DE*TEK618)  T I  ME ( N I -5 , UN -SE C )  P(UN- ' RAO l ANS/SEC ' , N l -3  ,  .  .  . 
SC«2  .  ,  LO- *  2 .  )  Q(LI»3,NI-3,SC-2.,LO»-2.)  R(LI-4,NIa3,  ... 
SC-2. , LO--2 . ) 

LABEL  G2  ROLL,  PITCH,  YAW  RESULTING  FROM  STEP  PITCH  INPUT- 

END 

STOP 


2.  Complete  Eaiations  of  Motion 

The  "following  source  code  is  the  DSL  program  used  to  generate  the 
results  discussed  in  Chapter  II  for  the  acnplete  AROD  model.  Notation 
follows  that  of  Appendix  ADA3A  wherever  possible.  Data  listed  in 
Appendix  ADATA  is  also  reproduced  here  in  the  form  used  toy  DSL. 


TITLE  AROD  EQUATIONS  OF  MOTION  -  FULL  MODEL,  NO  CONTROL  -  7  MAY  87 
CONST  GRAV-  32.174,  PI-  3.1415962 
ARRAY  T RANS ( 9  )  ,  ABTRAN ( 9 ) 


r  “  /  P  /  Y  D  E  L  « 
R/P/YSLOPE- 
THRMI N/MAX- 
R/P/JLEQ  « 
VEQ  > 

VANEFF  - 
CL/S/DDEL  ■ 
L/S/DSLOPE- 
FRCSlG  « 
WTRAT  .  - 

P/Q/R  « 

UU/VV/.WW  « 
PHI  /  THT/PSI 
ALPHA/BETA/A 
DELE/R/A  ■ 

THROT  ■ 

ALT  ■ 

U/V/WERTH  ■ 
GRNSPD/D I  ST 
THRHOV  - 
HOVRAT  ■ 

THR  a 


MNEMONICS  * 

ROLL/PITCH/YAU  COEFFICIENTS  * 

ROLL/PITCH/YAW  SLOPES  OF  COEFFICIENT  CURVES  * 

THRUST  FORCE  LIMITS  (LBS)  * 

EQUIVALENT  ROLL/PITCH/YAU  FORCE  FACTOR  * 

EQUIVALENT  VELOCITY  FACTOR  * 

CONTROL  SURFACE  EFFECTIVENESS  * 

L1FT/SIDEF0RCE/DRA6  COEFFICIENTS  * 

LI FT/SIDEFORCE/DRAG  SLOPES  OF  COEFFICIENT  CURVES  * 
CONVERSION  FOR  I N** 2 • L 8 ( FORCE )  TO  FT**2-SLUGS  * 

RATIO  OF  AROD  WEIGHT  TO  WIND  TUNNEL  MODEL  WEIGHT  * 

ROLL/PITCH/YAU  ANGLE  RATES  FOR  EOM  (RAD/SEC)  * 

BODY-FIXED  VELOCITY  COMPONENTS  (X/Y/Z )( FT/SEC )  * 

-  EARTH-FIXED  EULER  ANGLES  (RADIANS)  * 

OATOT--  ANGLES  OF  ATTACK  (RADIANS)  * 

ELEVATOR/RUDDER/AILERON/THROTTLE  POSITIONS  BASED  ON  * 
IDENTICAL  SERVO  EQUATIONS  (RADIANS)  * 

CHANGE  IN  THRUST  FORCE  BASED  ON  THROTTLE  SERVO  (LBS)* 
ALTITUDE  BASED  ON  INTEGRATION  OF  VERTICAL  SPEED  (FT)* 
EARTH  FIXED  VELOCITIES  (X/Y/Z)  (FT/SEC)  * 

-  SPEED  AND  DISTANCE  (NO  DIRECTION)  OVER  GROUND ( F T / S  )  * 

THRUST  FORCE  REQUIRED  IN  HOVER  (LBS)  * 

ANGULAR  RATE  OF  PROPELLER  IN  HOVER  (RAD/SEC)  BASED  * 
ON  EXPERIMENTAL  DATA  * 

TOTAL  THRUST  FORCE  (LBS)  * 
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NEUVEQ  ■■  I  HP  I R I  CAL  DATA  FOR  VELOCITY  CALCULATIONS  * 

VDELTA  ■  •  EQUIVALENT  VELOCITY  FOR  FORCE  CALCULATIONS  (FT/SEC)  * 

RHOVA  »•  WEIGHTING  FACTOR  FOR  AERO  COEFFICIENTS  BASED  ON  * 

VEQ  AND  AOATOT  * 

DELTIP  ■  »  CHANGE  IN  ROTOR  ANGULAR  VELOCITY  (RAD/SEC)  * 

VTIP  --  ROTOR  ANGULAR  VELOCITY  (DISCRETE)  (RAD/SEC)  * 

USEDEL/SLP/EQ--  COMBINATIONS  OF  VARIABLES  WITH  WHICH  AERODYNAMIC  * 
DATA  CAN  BE  USED  FOR  MOMENTS  AND  FORCES  (UNITS  BASED  ON  DATA  * 
L/S/DAF  ■■  LIFT/SIDE  FORCE/ORAG  AERO  FORCE  COMPONENT  (LBS)  * 

L/S/DTF  «■  LIFT/SIDE  FORCE/DRAG  THRUST  FORCE  COMPONENT  (LBS)  * 
FAX/Y/Z  =  «  AERODYNAMIC  FORCES  (LBS)  * 

FTX/Y/Z  --  THRUST  INDUCED  FORCES  (LBS)  * 

R/P/YAN  «•  ROLL/PITCH/YAW  AERODYNAMIC  MOMENT  COMPONENT  (FT-LBS)* 

R/P/YTM  --  ROLL/PITCH/YAW  THRUST  MOMENT  COMPONENT  (FT-LBS)  * 

LA/MA/NA  »■  ROLL/P  I TCH/YAW  AERODYNAMIC  MOMENT  (FT-LBS;  * 

LT/MT/NT  ■  ■  ROLL/PITCH/YAW  MOMENT  DUE  TO  THRUST/ROTOR  (FT-LBS)  * 

GX/Y/Z  *«  BODY-FIXED  GRAVITY  COMPONENT  (FT/SGC**2)  * 


EXPERIMENTALLY  OBTAINED  DATA 


UNITS 

* 

MOMENTS  OF  INERTIA 

I XX , I Y Y , I ZZ , I RX 

LB-FT-SEC**2 

FORCES 

FATX , FATY, FATZ.ETC 

LB 

MOMENTS 

LAT, MAT, NAT. ETC 

LB-FT 

ANGLES 

PSI  ,THT,PHI 

RADIANS 

PSIDEG, THTDEG, PHIDEG 

DEGREES 

ANGLE  RATES 

P.O.R 

RADIAN  / 

SEC 

* 

VELOCITIES 

UU, VV.WW.UERTH.ETC 

FT  /  SEC 

— 

UR 

RADIAN  / 

SEC 

CONST  IXXM-7063 .39, I  YYM-7768 . 22 , l ZZM-7729 . 58 , IRXM-69.552,  WEIGHT-  76.5 

*  MEASURED  IN**2-LB  IN**2.-LB  IN**2-LB  IN**2-LB  LB 

* 

CONST  PROPUT*  4.50017,  PROPEF-  10.472,  RHO-  .00192 

*  L  3  UNITLESS  PROPELLER  WE  I G H T/ E F F  I  C  I  E N C Y 

CONST  THRMIN-  35.,  THRMAX-  115. 

*  LB  LB  MIN/MAX  THROTTLE  POSITION 

CONST  TAMAX-  100.,  RMAXT *  100. 

*  LBS/SEC  LB/SEC*SEC  MAX  CHANGE/RATE  OF  CHANGE  OF  THRUST 
CONST  MAXDFL"  .5236,  RMAX-  .87266 

*  RADIANS  RAD/SEC  CONTROL  VANE  MAX  DEFLECTION  &  RATE 

CONST  LAEFF-  -21.29,  MEEFF-  -14.51,  NREFF-  -16.68 

*  ROLL,  PITCH,  YAW  VANE  EFFECTIVENESS 

CONST  HI-  17.77,  H2-  157.91,  HG1-  219.912,  HG2-  24674.126 
*ELEV,  RUD,  AIL,  THROTTLE  GAINS/  GYRO  GAINS 


CONST  WE-  1 


K  -  .5 


*  THROTTLE  INTEGRATION  LOOP  GAINS 

CONST  ALPN I NB  .8727,  ALPMAX*  1.570796,  AOAN I N«  50.,  AOANAX>  90. 

*  ANGLE  OF  ATTACK  LIMITS  FOR  TABLE  DATA  REFERENCE 


AERODYNAMIC  PARAMETERS  FROM  WIND  TUNNEL  TESTS 


AFGEN 

AFGEN 

FGEN2 


CYDEL 

YSLOPE 

YEP 


AFGEN 


VANEFF 


FGEN2 


VEO 


*  NOTE 

:  ALL 

TABLE  DATA 

IS  A  FUNCTION  OF  ANGLE 

OF  ATTACK 

(TOTAL) , 

* 

AOATOT 

RANGES  FROM  50  TO  90 

DEGREES . 

»*  * 

*  ROLL  DATA 

AFGEN 

CRDEL 

a 

50, .0807, 

55  ,  .0807, 

60, .0807, 

65, .0807, 

,  (  , 

70, .0743, 

75,  .071  2, 

80, .0810, 

90, .0850 

AFGEN 

RSLOPE 

a 

50, .1514, 

55, .1533, 

60, .1673, 

65, .1977, 

»  •  • 

70, .1834, 

75, .1764, 

80, .2041 , 

90, .2000 

FGEN2 

REG 

a 

8,2,16. 

50,55,60,65, 

70,75,80,90 

.  0,10,  . 

•  • 

9.293,8.5 

,  8.82,7.1, 

8. 115,7.2 

.  7.848,6 

»  4  i  •  •  • 

7.347,5.9 

,  7.134,5.4 

,  5.011 ,4.67,  0.0, 0.0 

****** 

r  *  * 

*  PITCH  DATA  ******** 

AFGEN 

CPOEL 

a 

50 , .31 06, 

55, .3106, 

60, .3106, 

65, .3106, 

•  •  • 

70 , .3081 , 

75  ,  .3305  , 

80  , .351  5  , 

90, .0000 

AFGEN 

PSLOPE 

a 

50, .5633, 

55  .  .51  23, 

60, .6974, 

65, .7612, 

•  •  • 

70, .7603, 

75 , .8189, 

80, .8853, 

90, .9000 

F  G  E  N  2 

PEG 

a 

8.2,16. 

50,55,60,65, 

70,75,80,90 

.  0.10,  . 

•  • 

34.138,11.69,  31.002,10.45,  29.576,9.72,  26.738,9.44,. 
24.545,6.78,  20.427,4.66,  1 5 . 1 68 , • 0 . 53 ,  0.0, -5.92 

**********  YAW  data  ********************************, 

60. . 0024, 

80, .0120, 

60,  .0215, 

80. . 0302, 


50, .0024, 
70, • .0015, 
50,  .0609, 


65. . 0024, 
90, .0000 

65. . 0060, 
90,  .0400- 

0.10,  . 
*  .  46 , 
*1.1, 


.  09 , * . 49 , . . . 

1 .8,  *1 .2 


55  ,  .0024, 

75,  .0090, 

55, .0663, 

70, -.0038,  75, .0223, 

8,2,16,  50,55.60,65,70,75,80,90, 

0.555,0.04,  0.576, -.51,  0.006, 

•1.634, -1.3,  -1.771, -.55,  -1.736, 

*  CONTROL  SURFACE  EFFECTIVENESS 

50.1.5,  55,1.5,  60.1.5,  6T.1.5 

70.1.5,  75,1.5,  80,1.9,  90,1.0 

***  VELOCITY  COMPENSATING  FACTOR 

8,2,16,  50,55,60,65,70,75,80,90, 

51.616,44.93,  44.955,39.11,  38.564,33.00,  32.935,27.45,. 
27.757,21.31,  21.829,15.06,  13.314,11.65,  0.0, 0.0 


***m 

0,10, 


Hitt 

LIFT  DATA 

AFGEN 

CLDEL 

m 

50,  .2428, 

55  ,  .  2428, 

60, .2428, 

65, .2428 

70, .2391 , 

75,  .1907, 

80, .1734, 

90, .1046 

AFGEN 

LSLOPE 

u 

50,  .7954, 

55,  .71  83, 

60, .6598, 

65, .5949 

70,  .5901 , 

75  ,  .4726, 

80, .4366, 

90, . 0000 

CONST 

LEQ 

a 

I  *  1 

85  . 

DRAG  DATA 

******** 

AFGEN 

CDOEL 

8 

50  ,  .5754, 

55  ,  .5754, 

60, .5754, 

65, .5754 

70,  .6052, 

75  ,  .5922, 

80, .4966, 

90, .3890 

AFGEN 

OSLQPE 

s 

50,1.3017, 

55,1.3131 

,  60,1.3711 

,  65,1. 

70,1,4933, 

75,1. 4674 

,  80,1.2507 

.  90,1. 

CONST 

DEG 

a 

0.0 
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**************************  SIDE  FORCE  DATA  **************** 
AFGEN  CSDEl  *  50, .0091,  55, .0091,  60, .0091,  65, .0091,  ... 

70, .0258,  75, .0331,  80, .0318,  90, .0000 

CONST  SSLOPE  *  0.0801 
CONST  SEQ  ■  0.0 


INITIAL 

* 

•**  CONVERSIONS 

*****  FORCE  TO  SLUGS  ANO  IN**2  TO  FT**2 
FRCSLG*  1 . / ( GR AV* 144.) 

IXX*  I XXM*FRCSLG 
IYY*  I YYM*FRCSLG 
IZZ*  I ZZN*  FRCSLG 
IRX*  IRXN*FRCSLG 
***  MOMENTS  OF  INERTIA 
IYMIZ*  IYY  •  IZZ 
IZNIX*  IZZ  -  IXX 
IXMIY*  IXX  ■  IYY 

* 

MASS*  WEIGHT /GRAV 
WTRAT*  WEIGHT/85. 

**********************  INITIAL  conditions  ************** 
PO  *  0 . 

QO  *  0 . 

RO  *  0. 

UO  *  0 . 

VO  *  0. 

WO  *  0. 

PHIO-  0. 

THTO*  0. 

PSJ.0-  0. 

DELEDO*  0. 

DELEO  -  0. 

DELRDO*  0. 

DELRO  *  0. 

DELADO*  0. 

*  DELAO  *  0. 

DLTRDO*  0. 

DLTRO  *  0. 

THROTO*  0. 

ALTO  «  0. 

GRNSPO*  0. 

THRHOV*  WEIGHT 

HOVRAT*  (THRHOV  -  PROPWT )*PROPEF 
************  DEFINE  MACRO  TO  LIMIT  ANGLE  *************** 
MACRO  ANGOUT*  AN L  I  M ( AN G I N , AN GM I N , AN GM AX  ) 

ANS I GN*  S  I  GN ( 1  .  ,  ANG  I  N  ) 

ANGOUT*  LIMIT(ANGMIN,ANGMAX,ANSIGN*ANG!N) 

ANGOUT*  ANSI GN*ANGOUT 

ENOMAC 

*********  DEFINE  MACRO  FOR  USE  IN  SERVO  EQUATION  LIMITING 
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MACRO  OUT  ■  SERVO(OUTO,OUTD,LIM) 
OUT"  IMTGRL(OUTO,OUTO> 

OUT"  LIMIT(*LIM,LIM,OUT) 

ENOMAC 


DYNAMIC 

* 


.***•**••*••  COORDINATE  TRANSLATION  COMPONENTS 
>«*  ROOT -FIXED  TO  EARTH  FIXED  (EULER  ANGLES) 
ROW  1 

TRANS(I)"  COSCPSI >*COS(THT) 

T  R£N  S  (  2  )  ■  COSCPSI  )*SIN(THT)*SIN(PHI  ) 

TRANS (3)-  COS (PS  I )*S I N ( T M T ) *COS ( PH I ) 

ROW  2 

TRANS(A)"  S I N ( PS  I  )*COS(THT) 

TRANS(5)>  S I N ( PS  I  )*SIN(THT)*SIN(PHI ) 

TRANS ( 6) ■  S I N ( P S  I  )*SIN(THT)*COS(PHI ) 

ROW  3 

TRANS(7)-  -SIN(THT) 

TRANS(8)«  COS(THT)*SIN(PHI )  . 

TRANS(9)>  COS(THT)*COS(PHI ) 


S I N(PS I )*COS(PN I ) 
SIN(PSt )*SIN(PHI) 


COSCPSI )*COS(PHI  ) 
SINCPHI )*COS(PSI ) 


*****************  VELOCITIES  *************** 
VTOT  •  SQ R T ( UU*UU  *  VV*VV  *  WU*WU) 

UERTH"  UU*TRANS( 1 )  ♦  VV*TRANS(2)  *  WW*TRANS(3) 
VERTH-  UU*  T  R AN  S ( 4  )  *  VV*TRANS(5)  *  WU*TRANS(6) 
UERTH"  UU*TRANS (7)  *  VV*TRANS(8)  ♦  UU*TRANS(9) 
GRNSPD"  SORT (VERTH*VERTH  *  UERTH*WERTH) 


********************  ANGLE  OF  ATTACK  **** 
I F ( VTOT . LE . .5}  THEN 

_  ALPHA"  1.570796 

BETA  >  0. 

ELSE 

ALPHA"  AS  I N ( UU/VTOT ) 

BETA  *  AS  I N ( VV/VTOT  ) 

ENOIF 

ALPHA  ■  AN L I N < ALPHA , ALPN I N , ALPNAX ) 

BETA  •  ANLIMCBETA.O. , 1 .570796) 

AOATOT"  RADEG*SQRT (ALPHA*ALPHA  *  BETA* BETA) 
AO A TOT"  LIMIT(AOAMIN, AO AM AX , AOATOT  ) 


,***********  ANGLE  OF  ATTACK  COORDINATE  TRANSFORMATION 
>  BODY-FIXED  TO  VELOCITY 
ROW  1 

ABTRAN ( 1 )■  -COS(ALPHA)*COS(BETA) 

ABTRAN ( 2 ) ■  SIN(BETA) 

ABTRAN (3  )■  S I N ( A L PH A ) *COS ( BE T A > 

ROU  2 

A8TRAN(4)«  COS ( A L P HA ) *$ l N ( BE T A ) 

ABTRAN ( 5 ) >  COS(BETA) 

ABTRAN(6)«  • S  I  N ( ALPH A ) *S I N ( BE TA ) 
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*  ROW  3 

ABTRAN ( 7) ■  • S 1 N ( ALPHA ) 

ABTRAN(8)>  0. 

ABTRAN<9)>  -COS(ALPHA) 

* 

••****••»*••••**»**•  FORCE  components  ***************************** 
THR  ■  L  I  M  I  T ( T  HRM I N  ,  THRNAX ,  THRHOV*THROT ) 

NEWVEQa  FGEN2(VEO, AOATOT, 0.  ) 

VDELTA*  VTOT*  SORT (WTRAT)*NEWVEO 
RHOVA  «  RHO*P!*SORT (UTRAT )*NEWVEO 

* 

DE LT  I  P»  (THR  *  WEIGHT)*PROPEF  _ 

VtIp  »  HOVRAT  ♦  OELTIP 

* 

USEDEL"  •DELTIP*RHOVA*UTRAT 
USESLPa  VDELTA*UTRAT 

USEEQ  ■  Cl.  ♦  (2./VTIP)*DELTIP)  *  WTRAT 
***  FORCES 

LAF  >  USEDEL*AFGEN(CLDEL, AOATOT)  ♦  USE SLP*A F GEN ( L S LOPE , AOATOT ) 
SAF  «  USEDEL*AFGEH<CSOEL, AOATOT)  ♦  USESLP*SSLOPE 
OAF  >  USEDEL*AFGEN(CDDEL, AOATOT)  ♦  USE S L P* A F GE N < D S L OP E , AO AT OT ) 
LTF  «  USEEQUEQ 

* 

FAX  •  ABTRANC 1 )*DAF  ♦  ABTRAN (2  )  *S A F  ♦  ABT R AN ( 3 ) * L A F 

FAY  •  ABTRAN (4  )*OAF  ♦  ABT R AN ( 5  ) * S A F  ♦  ABT R AN < 6  ) * L A F 

FAZ  a  ABTRAN(7)*DAF  ♦  AB T R AN < 8 ) * S A F  ♦  AB T R AN ( 9 ) * L A F 

FTX  a  AB T R AN ( 3  )  * L T F 

FTY  ■  ABTRAN(6)*LTF 

FTZ  •  ABTRAN(9)*LTF 

FATX  a  FAX  ♦  FTX  . 

FATY  a  FAY  ♦  FTY 

FJLtZ  ■  FAZ  ♦  FTZ 

***  MOMENTS 

RAM  a  USEDEL*AFGEN(CROEL, AOATOT)  ♦  USESL P*A FGEN < RS LOPE , AOATOT )  . 
+  USEEO*FGEN2<REQ, AOATOT, 0. ) 

PAM  a  USEOEL*AFGEN<CPOEL, AOATOT)  ♦  US E S LP* A FGEN ( PS  LOPE , AOATOT )  . 
+  USEEQ*FGEN2(PEQ, AOATOT, 0. ) 

YAM  a  USEDEL*AFGEN<CYDEL, AOATOT)  +  USE S L P*A F GEN ( YS LOPE , AOATOT )  . 
♦  USEE0*FGEN2(YE0, AOATOT, 0. ) 

* 

LA  •  ABTRAN( 1 )*RAN  +ABTRAN ( 2 ) *P AN  ♦ABTRAN ( 3 )*YAM 
MA  >  ABTRAN(4)*RAM  ♦ABTRAN ( 5 )*PAH  ♦ABTR AN ( 6 ) * Y AN 
NA  a  ABTRAN(7)*RAM  ♦ABTR AN < 8 )*PAN  ♦ ABT R A N ( 9 ) * Y AM 
LT  ■  -DELT1P*DELTIP*FRCSLG*.0729  ♦  I  XX*LAE F F*DELA 
NT  «  MEEFFM  YY*DELE 
NT  >  NRE  F  F*  1  ZZ*DELR 

LAT  ■  LA  ♦  LT 

MAT  ■  MA  +  NT 

NAT  »  NA  ♦  NT 
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GRAVITY  COMPONENTS 


GX  -  -GRAVMRAHSd  > 

GY  -  -GRAV*< *TRANS<4)) 
GZ  -  •GRAV*(*TRAMS<7)> 


SAMPLE 

***»  INSERT  CONTROLLER  HERE 


DERIVATIVE 

* 

a*******************  MOMENT  derivatives  ******** 
PD  ■  ( I YM I Z*R*Q  •  1 RX*DELT I P  ♦  LAT)/IXX 
QD  *  < I ZM I X*P*R  •  I RX* VT I P*R  ♦  MAT  )  / I Y Y 
RD  *  ( I XM I Y*P*G  ♦  I RX*VT I P*Q  ♦  NAT)/IZZ 
*•*******•***•*••*••  FORCE  derivatives  ********* 
UD  •  (VV*R  •  UU*Q  >  *  FATX/MASS  *  GX 
VD  -  <VU*P  -  UU*R )  *  FATY/MASS  ♦  GY 
UD  *  ( UU*Q  -  VV*P)  *  FATZ/MASS  ♦  GZ 

•••••••••A**********  ANGLE  derivatives  ********* 

PNID-  P  *  ( 0*S I N ( PH  I  )  *  R*COS(PHI  ))*TAN(THT) 
THTD*  Q*COS(PHI)  •  R*S I N ( PH  I ) 

PSID-  (QMIN(PNI)  *  R*COS(PHI  ))/COS(THT) 
*••••«••••****•*••••  SERVO  DERIVATIVES  ********* 
***  ELEVATOR,  RUDDER,  AILERON,  THROTTLE 
DELEDD*  •  H 1 *D  E  LED  •  H2*DELE  *  H2*UE 
DELRDD*  -  H  1  *D E LRD  •  H2*DELR  *  H2*UR 
DELADD*  *.H  1  *DELAD  •  H2*DEL A  *  H2*UA 
OLTRDD*  •  H1*DLTRO  *  N2*DLTR  ♦  N2*UT 
tHJLOTD  »  •  UE*T  HROT  ♦  K*UE*DLTR 


*** INTEGRATIONS 


••••••*••••••**•••  MOMENT  EQUATIONS 

P  *  INTGRL(PO.PD) 

Q  *  INTGRL(QO.QD) 

R  *  INTGRL(RO.RD) 

*•*•««•••••••••*•••  FORCE  equations 

UU ■  INTGRL(UO.UD) 

Wa  I  NTGRL  (  VO  ,  VD  ) 

UU-  INTGRL(UO,UD) 

a*******************  HULER  angles 
PHI  a  INTGRL(PHIO.PHID) 

THT  -  INTGRL(THTO,THTD) 

PSI  -  INTGRL(PSIO.PSID) 

PHIDEG-  RA0EG*PH  I 
THTDEGa  RADEG*THT 
PSIDEG-  RADEG*PS! 


a*******************  SERVO  EQUATIONS 
DELED  >  SERVO(DELED0,DELEDD,RMAX) 
DELRD  -  SERVO(DELRDO, DELRDD, RMAX) 


a 
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DELAD  *  SERVO(DELADO,DELADD,RNAX) 

DLTRD  ■  SERVO(DLTRD0,DLTRDD,RNAXT) 

DELE  ■  INTGRUDELEO, DELED) 

DELR  ■  SERVO<DELRO,DELRD,NAXDFL) 

DELA  ■  I  NTGRLCDELAO, DELAD) 

DLTR  ■  SERVO(DLTRO, DLTRD, TANAX) 

THROT  ■  I HTGRL(THROTO, THROTD  ) 

***  ALTITUDE 

ALT  ■  I NTGRL CALTO , UERTH  ) 

D I  ST*  1MTGRL(GRNSP0,GRNSPD) 

CONTRL  F1NTIN-  5 

PRINT  .5^  PHI  DEG, THTDEG.PS 1  DEG, P,Q,R,UU,VV,WW,DELA 
SAVE  .05,  PHIDEG,  T  HTDEG  ,  PS  I  D  E  G  ,  P  ,  Q  ,  R  ,  II U  ,  VV  ,  VW 
PARAN  UE  *0 .  , UR  *0 . , UA- • . 052 , UT *  O.,DELAO*0 
LABEL  ( G 1 )  TINE  RESPONSE  TO  UA  --.052(NO  CONTROL) 

GRAPH  ( G 1 ,  DE-TE  K61. 8  )  T I  ME ( N I *5  ,  UN ■ SEC  )  P H I DE G ( N I *4 , SC-45 , L I  *  1 
UM-'THTDEG  AND  PSIDEG  ( DE GREE S  )  ' , LO« • 90 )  T HTD E G ( AX-ON  I T , N I *4 
SC *45, LI-3, LO>- 90)  PSIDEG(AX- ON  IT, NI-4, SC-45, LI-4, L0*'90) 
GRAPH  (G2,DE*TEK618,0V,PO»0,S)  T I NE ( N I -5 , UH-SE C  )  P(NI-4,  ... 

UH-'O  AND  R  (RAD/SEC)  •  ,LO--2, SC-1 , LI-1 )  Q(AX-ONIT,  ... 
HI -4, LO-- 2.SC-1 ,LI»3)  R ( AX-ON  I T , N I -4 , LO- -  2 , SC- 1 , L I -4 ) 
GRAPH  <G2,DE-TEK618,0V,P0-6,Q)  T I NE (N I -5 , UN-SEC )  VV(NI-4,  ... 
UN- 1 UU  AND  UU  (FT/SEC)' ,LO-*2, SC-5,  LI-3)  UU(AX-0NIT,  .. 
Nl-4,LO-*2, SC-5, LI-1 )  H«< AX-ON  I T , N I -4 , LO- -  2 . SC-b , L I -4 ) 

END 

STOP 
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AHSMULF  C 

axnsauai  aid  linear  keel  snnuenoi 

Hie  DSI/VS  source  code  used  to  simulate  the  linearized  model  is 
listed  here.  Hie  steady-state  controller  is  inplenented  under  the  ~ 
"SAMPLE"  segment.  This  same  controller  is  used  in  the  nonlinear  model 


TITLE  SIMPLE  MODEL  WITH  OPTIMAL  GAINS  FROM  CONTROLS  -  9  AUG  87 

***  FEED  BACK  GAINS  USING  THE  ERROR  STATE  CONTROL  LAW  U«  -F(X*R)  *** 

ARRAY  R  KPT ( 8 ) , RKQR ( 8 ) ,  XP T (8) ,708(8) , UP  T ( 2 ) , UQR ( 2 ) , EPT ( 8) , EQR ( 8 ) 

FGEN2  FBKPT «  2,8,16,  1 , 2 , 1 , 2 , 3 , A  ,  5 , 6 , 7 , 8  ,  ... 

•  1.99484,- .60299. - . 71 71 7 , . 27437 , 1 . 74559 , . 00994 , . 1 4049 , . 000 1 6 ,  . . 
-.58535,1  1  1  .  16428. • .20533,28.90878, . 581 70 , 1 . 66845 ,. 00962 ,.  1 3667 

FGEN2  FBKQR-  2,8,16,  1 , 2 , 1 , 2 , 3 , 4 , 5 , 6 , 7, 8 ,  ... 

•1  .61288, -  1 .61027, - .55320, .31985, 1 .05412, • .20752, .  12656, - .00264,  . 
1.56829,- 1.59987, • .30479, • .56615. . 1 0693 , 1 . 23477 , . 00096 , .12962 

* 

CONST  RMAXT ■  100.,  TANAX-  100.,  NAX0PL-  .5236,  RNAX-  .87266 
CONST  Ala. 0236133, A2*-3. 73 81  , A3» * 6 . 750686 , A4a . 4830933 
CONST  A5j»6. 784433, A6«. 107891  ,A7». 893727 
CONST  LAE«- 21 .29,NEE«- 14.51 .NRE-* 16.68 
CONST  H 1 ■  17.77,  H2a  157.91 

CONST  PHDES-.0,ALRDESa.00,THDESa.0,PSDESa.0 

CONST  RKPT<3)a0,RKPT<4)a0,RKPT(5)a0,RKPT<6)a0,RKPTC7)»0,RKPT(8)a0 
CONST  RKQR(3)aO,RKQR(4)aO,RKQR(5)aO,RKQR(6)aO,RKQR(7)-0,RKQR(8)aO 
INIT 

DELEDOa  0 
DELRDOa  0 
D E LAD 0a  0 
DLTRDOa  0 
DELEOa  0 
DELR0-  0 
DELA0-  0 
OLTRQa  0 
P0a  o 
Q  0  ■  0 
ROa  0 
THROa  o 

I  a  1 
J  ■  1 
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*********  DIMM  MACRO  FOR  UII  IN  tIRVO  EQUATION  LIMITING 

MACRO  OUT*  SRRVO<OUTO,OUTO,Ll.M) 

OUT*  INTGRL(OUTQ,OUTO) 

OUT*  LINIT<-llN,LIM,OUT> 

CNOMAC 

DYNAMIC 

****  INRUT  CALCULATIONS 

RKRT(I)  *  PHDES*STER(0.  ) 

RKRT(2)  *  ALRDES*STSR(Q. ) 

RKOR(1)  ■  T  HD  E  S*S  T  ER  (  0  .  ) 

RKQR ( 2  )  *  RSDES*STER( 0 . ) 

*  ,  „ 

XRT( 1 >■  RNIS 
XRT ( 2 ) »  U 
XRT (3 ) ■  R 
XRT ( A ) *  THROT 
XRT  <  5 ) *  DELA 
XRT ( 6 ) *  DLTR 
XRT (7>*  DELAD 
XPTC8)*  DLTRO 
XQR( 1  )■  THTN 
XQR ( 2 ) ■  RS1B 
XQR ( 3 ) *  Q 
XQR(A)*  R 
XQR ( 5 ) ■  DELE 
XQR (6 ) ■  DELR 
XQR ( 7  > ■  DELED 
XQR(8}«  DELRD 

SAMPLE 


DO — 30  1*1,2 
•URT  -0 
BUQR  *0 
DO  20  J«1 ,8 

ERT  (  J  )■  XRT  <  <1  >  -  RKRT  (  J  ) 
EQR ( J ) *  XQR ( J )  •  RKQR(J) 
FBPT*  F  GEN2 ( FBKRT , I , J ) 
FBQR*  F  GEN2 ( FBKQR , I , J ) 
BURT*  BURT  ♦  EPT(J)*FIPT 
BUQR*  BUQR  ♦  EQR(J)*FBQR 
20  CONTINUE 

URTU)*  -(BURT) 

UQR< I )■  -(BUQR) 

30  CONTINUE 
* 

UA*  URT ( 1 ) 

UT*  URT ( 2 ) 

UE*  UQR  (  1  ) 

UR*  UQR ( 2 ) 
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DER1V 


MUD*  P 
UD  ■  A7*TNR0T 
THTBD-  0 
MUD*  It 

PD-  A 1  *U  ♦  A2*THR0T  ♦  LA(*0«L A 
T  HROTD*  • THROT  ♦  .  5 * 0 L T R 
QD»  A3*R  ♦  NRK*DELR 
RD*  A5*8  ♦  NRE*DBLR 

*»*  ILRVATOR ,  RUDOIR ,  A1LRR0R,  THROTTLE 
DELADD-  • H1*DELAD  •  H2*DELA  ♦  N2*UA 

DLJRDD-  *H1*DLTRD  •  H2*0LTR  ♦  H2*UT  _ 

DKLEDDa  -N1*DELED  •  N2*DELE  ♦  N2*UE 
DELRDDa  >H1*DELRD  *  N2*DELR  ♦  H2*UR 

* 

PHIB-  INTGRL(PHIBO.PHIBD) 

U  ■  1  NTGRL (UO , UD } 

THTIa  IHTGRL(THTBQ, THTBD) 

Mila  INTGAL(PS t BO.PSIBD) 

Pa  INTGRL(PO.PD) 

THROT-  INTGRL(THRO.THROTD) 

Q-  INTGRL(QO.QD) 

Ra  INTGRL(RO.RD) 

* 

i  DELED  «  SERVO(DELEDO,DELEDD,RNAX) 

j  DELRD  >  SERVO(DELRDO,DELRDD,RNAX) 

DELAD  •  SERVO(DELADO,DELADD,RNAX) 

DLTRD  -  SERVO(DLTRDO,DLTRDD.RNAXT) 

DELE  a  INTGRLCDELEO, DELED) 

DELR  >  $ERVO(DELR0, DELRD, HAXDPL) 

DELA  -  IHTGRUDELAO,  DELAD) 

a  SERVO(DLTRO,DLTRD,TANAX) 

DESa. 1 
ZERa  o. 

CONTROL  FlNTlN-5,DELS-.04,DElT-.005 
PARAN  PHIBO-.OO.UO-Q.TNTBO-O.PSIBO-O 
SAVE  .04,PHIB,U,TNT3,PSIB,ZER,DES 
l  PRINT  .  2,PHII,U,THT1,PSU 

»  LAIEL  (61)  REGULATOR  GAINS  -  ERROR  STATE  •  PHI  (0)a  .1 

(  GRAPH  (61 ,DEaTEK618)  T I NE ( N I «5 , UN-SEC )  PH  I B ( N I -4 , SC- . 05 , L ! ■ 1 .  ... 

UN-'R/SEC  AND  ALT  CATE  F/ SEC • , LO- * . 05 )  U ( AX-ON  I T , N ! -A ,  ... 

SC-.S, LI-3, LO-- .5)  DES(AX-ONlT,Nl-10,SC-.05,U-4,lO-- .05)  ... 
ZER(AX*ONIT,NI«4,SCa. 05, LI-4, LO-* .05) 

GRAPH  (62,DEaTEK616,OV,POa6,4.5)  T  I  ME  (  N I -5  ,  UN-SEC  )  THTKNI-4,  .  .  . 

UN-'AND  PSIB  R/SEC' ,LO»- .05, SC-. 05,11-1)  PS  I  B ( AX-OH I T ,  ... 
NI-4.L0--  .05,SC>.05,LI>3)  0 E S ( AX-OH  I  T , N I -4 , LO- • . 05 , SC- . 05 , L I -4 )  . . . 
*  ZER(AX»OHIT, HI-4, SC-. 05, LI-4.L0-* .05) 

l  *ND 

t  STOP 
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H  ,  — to 

Tha  Numerical  Elistranagnati.es  Code  (NEC)  enables  the  uaar  to 
simulate  a  solid  surfaoa  using  either  tha  wira  grid  or  tha  surfaoa  patch 
technique.  Extensive  use  of  symmetry  in  the  program  and  the  circular 
symmetry  of  the  ARDD  model  allow  that  the  surfaoa  patch  approach  be 
easily  implemented.  Figure  D.l  illustrates  the  procedure  schematically. 


Figure  0.1  Schematic  for  ARDD  Body  Modelling 
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The  following  ocwputer  coda  is  the  NBC  nodal  of  the  body  which  was 
utiliaad. 

(1)  Build  the  1st  oolum  of  patches.  The  column  is  9.92°  in  angular 
width  (arc)  and  14  inches  high.  The  arc  of  the  oolunn  is  in  the 
xy  plana  and  the  height  of  the  oolunn  is  in  the  Z  direction  at  a 
radius  of  15  inches  from  the  Z  axis. 


SR 

0. 

o. 

IS., 

0.,  S..  18.435, 

o.. 

9.92 

SR 

0. 

o. 

15.. 

0.,  1.5,  5.711, 

o.. 

7.44 

SR 

0. 

o. 

15., 

0., -1.5, -5. 711, 

o.. 

7.44 

SR* 

0. 

o. 

15.. 

0.,- 5.  ,*18. 435, 

o.. 

9.92 

(2)  Rotate  and  duplicate  the  column.  The  1st  oolunn  is  rotated  in  the 
X?  plane  and  d^licatad  to  fora  a  oolunn  that  is  now  18.84°  wide 
and  14  inches  high. 

GM  0#  1 1  0 • i  0 • )  9 .4735  #  0 . (  0 • #  0* 

(3)  Add  patches  to  top  and  bottom  to  form  wedge.  Thts*  patches 
together  with  the  oolunn  form  a  18.84°  section  of  the  cylinder. 


SR 

0. 

0. 

13.95, 

1.16, 

7., 

90.. 

4.737, 

9.26 

SR 

0. 

0. 

13.95, 

1.16, 

*7., 

•90.. 

4.737, 

9.26 

SR 

0. 

0. 

11.94, 

.990, 

7.. 

90., 

4.737, 

7.94 

SR 

0. 

0. 

11.96, 

.990, 

*7., 

*90., 

4.737, 

7.94 

SR 

0. 

0. 

9.470, 

.785, 

7.. 

90., 

4.737, 

9.425 

SR 

0. 

0. 

9.470, 

.785, 

*7., 

•90.. 

4.737, 

9.425 

SR 

0. 

0. 

4.980, 

.413, 

7.. 

90., 

4.737, 

9.92 

SR 

0. 

0. 

4.980. 

.413. 

-7.. 

•90., 

4.737, 

9.92 

(4)  Oocqplete  the  cylinder.  The  entire  18.84°  wedge  is  rotated  19 
tines  to  form  a  360°  cylinder  which  is  missing  only  a  top  and 
bottom  center  piece.  These  two  canter  pieces  are  added  along  with 
each  antenna  that  is  to  be  simulated,  gs  scales  the  dimensions  to 
meters  from  inches  and  ge  declares  an  end  to  geometry. 

GR  0.  19 
GS  2 
GE 

(5)  Specify  frequency.  The  frequency  at  which  the  model  and  antenna 
are  to  simulated  is  specified,  m  this  exasple,  f  -  800  Jffe. 

PR  o, 0,0,0, 800.  ,o. 

(6)  Store  the  model.  The  entire  symmetric  modal  is  stored  in  a 
Green's  function  matrix  for  recall  when  the  antenna  is  simulated. 

WG 
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A  dipole  array  and  a  loop  antenna  are  simulated  in  the  NEC  oode 
which  follows, 
a.  5  Dipole  Ring  Array 

A  five  dipole  ring  array  consists  of  five  half-wave  dipoles  that 
are  arranged  in  a  circle. 

(1)  Rood  the  (Steen's  function.  The  APOD  body  model  which  is  stored  in 
a  Green's  function  matrix  is  recalled  to  start  the  simulation. 

This  matrix  also  includes  the  frequency  at  which  the  antenna  will 
transmit. 

GF 

(2)  Specify  the  dipoles.  Wire  .1  inch  in  radius  and  7.38  inches  long 
is  used  to  model  a  half-wave  dipole  at  f  -  800  MHz.  11  segments 
in  each  dipole  is  used  so  that  NEC  automatically  calculates 
current  losses.  5  dipoles  are  spaced  equally  around  the  circle  at 
a  radius  of  15.5  inches.  The  dipoles  longitudinal  axis  is  in  the 
Z  direction. 

GU  1,  11,  15.5,  0.,  3.69,  15.5,  0..-3.69,  .1 
GW  2,  11,  4.79,  14.74,  3.69,  4.79,  14.74,-3.69,  .1 

GU  3,  11,  -12. S4.  9.11  ,  3.69,  -12.54,  9.11  ,-3.69,  .1 
GU  4,  11,  -12.54,  -9.11,  3.69,  -12.54,  -9.11,-3.69,  .1 
GU  5.  11,  4.79, -14.74,  3.69,  4.79,  -14.74,-3.69,  .1 

(3)  Complete  the  cylinder.  The  two  center  patches  which  were  excluded 
from  the  symmetric  AROD  body  model  are  added,  the  dimensions  are 
scaled  to  meters,  and  the  geometry  is  ended  with  no  ground  plane. 

SP  0,  0,  0. ,  0. ,  7, ,  90,  0. ,  12.57 
$P  0,  0,  0.,  0  . , -  7 . , -  90 ,  0.,  12.57 
GS  2 
GE  0 

(4)  Specify  the  driving  source.  In  this  example,  a  feed  point  at  the 
center  of  each  dipole  provides  1.0  volts. 

ex  o,  i ,  6,  01 ,  i . ,  o. 


b.  Loop  Antenna 

A  single  wire  loop  antenna  is  added  to  the  APOD  body  model.  The 
loop  is  flush  to  the  surface  of  the  APOD. 
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(1)  Read  the  Green's  function.  The  ABOD  bod/  model  which,  is  stared  in 
a  Green's  function  matrix  is  recalled  to  start  the  simulation. 

This  matrix  also  includes  the  frequency  at  which  the  antenna  will 
transmit. 


G  f 

(2)  Specify  1st  wire  in  the  circle.  The  circle  will  consist  of  wire 
arc  segments  connected  together.  The  1st  segment  spans  10°  and  is 
at  a  radius  of  15.5  inches.  The  wire  is  in  the  XZ  plane  initially 
and  is  then  moved  to  the  XV  plane. 

Ga‘i ,3, 15.5, -5., 5., .1 

GN  0,  0,  90.,  0.,  0.,  0.,  0.,  0.,  1 

(3)  Create  the  loop.  The  first  wire  is  duplicated  and  rotated  36 
times  to  form  a  ccnplete  circle. 

GR  1,  36 

(3)  Ccnplete  the  cylinder.  The  two  center  patches  which  were  excluded 
from  the  symmetric  ABOD  body  model  are  added,  the  dimensions  are 
scaled  to  meters,  and  the  geometry  is  ended  with  no  ground  plane. 

SP  0,  0,  0.,  0.,  7.,  90,  0.,  12.57 
SP  0,  0,  0. ,  0. ,  7. , -90,  0.  ,  12.57 
GS  2 
GE  0 

(4)  Specify  the  driving  source.  In  this  example,  a  feed  point  at  <p  = 
0°  provides  1.0  volts. 

EX  0,  1,  6,  01,  1.,  0. 
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